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Introduction 


Chemotherapy  has  been  used  with  other  cancer  therapies  to  increase  the  survival  rate  of  patients,  but  it  also 
causes  serious  side  effects  due  to  its  great  toxicity  for  tissues  with  high  growth  fractions  such  as  the  bone 
marrow,  gastrointestinal  epithelium,  hair  follicles,  and  gonadal  tissue.  For  the  breast  cancer  therapy, 
chemotherapy  has  been  used  as  an  adjuvant  therapy  or  a  neoadjuvant  therapy  (or  induction  chemotherapy). 
In  both  of  the  cases,  a  chemotherapeutic  agent  is  given  by  systemic  administration,  and  thus  side  effects  are 
often  occurred  to  breast  cancer  patients.  Therefore,  delivery  of  chemotherapeutic  drugs  only  into  breast 
cancer  cells  will  reduce  side  effects  of  chemotherapy  significantly,  and  the  patients  will  be  more 
comfortable  during  chemotherapy. 

Currently,  MRI  and  CT  are  the  only  tools  for  monitoring  the  effects  of  chemotherapy  by 
measuring  the  tumor  size.  However,  detecting  early  effects  of  therapy  is  needed  to  improve  therapy 
efficiency  and  to  increase  patients’  survival  rate.  In  the  last  decade,  significant  investigations  in  both 
laboratory  and  clinical  settings  using  near-infrared  spectroscopy  (NIRS)  have  been  conducted  for  non- 
invasive,  quantitative  measurements  and  imaging  of  tissue  oxygenation  of  the  brain  '  •  and  muscle  ’ ' 
in  vivo.  Recently,  Liu  et  al.<>  and  other  groups10 11  have  shown  that  the  NIR  techniques  can  be  also  used 
for  monitoring  vascular  oxygen  dynamics  of  breast  tumors  during  respiratory  intervention  non-invasively. 

We  have  previously  studied  tumor  oxygenation  under  gas  intervention  using  NIRS  and  found  that 
oxyhemoglobin  concentration  changes  during  gas  intervention  can  be  fitted  by  a  two-exponential 
equation  containing  two  time  constants.Error!  Bookmark  1,01  der,ned-  \ye  suggested  that  those  two  time  constants 
are  related  to  blood  perfusion  in  tumors.  If  a  drug  is  delivered  into  a  breast  tumor  by  drug-loaded, 
biodegradable  fibers,  the  drug  will  be  released  most  likely  into  the  tumor  cells,  and  this  will  decrease  the 
side  effects  of  chemotherapy.  Our  hypothesis  for  this  proposed  work  is  that:  when  the  drug  is  delivered 
locally  by  biodegradable  fibers  to  the  tumor  cells,  changes  in  blood  perfusion  in  the  tumor  will  occur  and 
can  be  observed  by  the  NIRS  measurement.  If  the  tumor  vascular  and  tissue  oxygenation  changes  during 
respiratory  challenges  show  different  and  significant  responses  before  and  after  the  administration  of  the 
chemotherapeutic  agent,  it  may  indicate  that  NIRS  could  be  a  novel  technical  means  in  monitoring  early 
effects  of  chemotherapy,  which  will  help  physicians  make  proper  tumor  treatment  plans  for  individual 
patients. 

Therefore,  this  study  will  allow  us  to  validate  a  non-invasive  tool  that  can  monitor  the  effects  of 
chemotherapy  on  breast  tumors  and  to  demonstrate  advantages  of  local  drug  delivery  using  biodegradable 
fibers  for  breast  tumor  treatment.  Our  study  may  help  breast  cancer  patients  under  chemotherapy  have  a 
higher  survival  rate  and  better  quality  of  life  with  reduced  side  effects  and  proper  tumor  treatment 
selection. 


Objectives: 

The  goals  of  this  research  are  to  evaluate  NIRS  technique  and  |l)F  MR  p02  mapping  as  non-invasive  tools 
for  monitoring  the  drug  effects  in  breast  tumors  and  to  compare  the  effects  of  systemically  delivered  drug 
and  locally  delivered  drug  with  embedded  biodegradable  fibers.  The  project  has  five  specific  aims  and 
tasks  as  shown  below: 

Aim  1:  to  monitor  tumor  vascular  oxygenation  during  oxygen/carbogen  interventions  using  both  NIRS 
and  l9F  MR  p02  mapping  after  a  systemic  drug  delivery. 
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Aim  2:  to  compare  the  effects  of  the  systemically  delivered  drugs  on  the  vascular  and  tissue  oxygenation 
of  the  breast  tumors  in  order  to  validate  the  monitoring  capability  of  the  two  techniques  for  breast  cancer 
treatment. 

Aim  3:  to  develop/implement  biodegradable  fibers  loaded  with  a  chemotherapeutic  agent  for  local  drug 
delivery  at  the  specific  breast  cancer  site. 

Aim  4:  to  monitor  tumor  vascular  oxygenation  during  oxygen/carbogen  interventions  using  both  NIRS 
and  19F  MR  p02  mapping  after  a  local  drug  delivery  through  a  biodegradable  fiber. 

Aim  5:  to  compare  the  chemotherapy  effects  between  local  drug  delivery  and  systemic  drug  delivery. 

Task  /:  To  obtain  the  knowledge  and  appropriate  training  for  handling  laboratory  animals  and  performing 
the  measurements,  (months  1-7): 

a.  Get  training  in  handling  laboratory  animals  (month  1). 

b.  Study  the  fundamentals  and  operation  of  the  NIR  oximeter  (months  2-3). 

c.  Study  the  fundamentals  and  operation  of  the  l9F  MR  EPI  of  HFB  method  (months  4-6). 

d.  Gain  basic  understanding  of  tumor  physiology,  oxygenation  and  dynamic  processes  (month 
7). 

Task  2:  To  monitor  tumor  vascular  oxygenation  during  oxygen/carbogen  interventions  using  both  NIRS 
and  19F  MR  p02  mapping  after  a  systemic  drug  delivery,  (months  8-16): 

a.  To  implant  breast  tumors  on  Fisher  rats  (month  8). 

b.  To  perform  NIRS  and  19F  MR  p02  mapping  measurement  on  breast  tumors  with 
oxygen/carbogen  inhalation  after  systemic  drug  administration,  and  evaluate  oxyhemoglobin 
level  and  tumor  tissue  p02  changes  from  NIRS  and  l9F  MR  p02  mapping,  respectively, 
(months  9-14) 

c.  To  measure  tumor  size  and  hemoglobin  level  in  blood  periodically  after  systemic  drug 
delivery  to  observe  the  effects  of  systemic  chemotherapy,  (months  9-14) 

d.  To  compare  the  effects  of  the  systemically  delivered  drugs  on  the  vascular  and  tissue 
oxygenation  of  the  breast  tumors  in  order  to  validate  the  monitoring  capability  of  the 
two  techniques  for  breast  cancer  treatment,  (month  15) 

e.  To  study  the  effects  of  the  systemically  delivered  drug,  e.g.  tumor  size  reduction  and 
time  constant  during  respiratory  changes,  (month  16) 

Body  of  the  Report 

As  mentioned  above,  the  purposes  of  this  project  are  to  deliver  a  therapeutic  agent  into  the  breast  tumor 
directly  by  biodegradable  fibers  and  to  utilize  NIRS  and  l9F  MR  p02  mapping  as  tools  to  monitoring  the 
therapeutic  effects  of  the  chemotherapeutic  agent.  In  the  first  year,  my  focus  was  to  obtain  the  appropriate 
knowledge  and  skills  for  conducting  the  project  (months  1-7),  and  to  perform  the  measurements  of 
vascular  oxygen  dynamics  in  breast  tumors  by  using  NIRS  system  during  systemic  cancer  drug  delivery.  I 
have  mainly  accomplished  these  tasks,  as  reported  below: 

Accomplishments  during  the  first  period  (months  1-7): 

a.  Obtained  the  skills  for  handing  laboratory  animals.  I  have  taken  the  animal  handling  class 
provided  by  the  University  of  Texas  Southwestern  Medical  Center  at  Dallas,  also  taken  the  online 
training  from  the  University  of  Texas  at  Arlington,  and  practiced  saline  injection  through  i.p.  and 
i.v. 

b.  Studied  the  fundamentals  and  operation  of  the  NIR  oximeter  (see  the  reportable  outcomes). 
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c.  Studied  the  fundamentals  and  operation  of  the  l9F  MR  EPI  of  HFB  method.  However,  this 
method  has  not  been  used  in  the  first  year  of  project  due  to  mismatch  of  the  animal  and  tumor 
implantation  schedule. 

d.  Gained  basic  understanding  of  cell  biology  and  tumor  biology.  I  have  taken  a  course  of  Cell 
Physiology  at  the  University  of  Texas  Arlington  and  also  Molecular  and  Cell  Biology  of  Cancer 
provided  by  the  University  of  North  Texas  Health  Science  Center  at  Fort  Worth  (spring  2004).  In 
addition,  I  have  taken  a  “Drug  Delivery”  class  at  the  University  of  Texas  Arlington  (fall  2003) 
since  I  will  use  the  local  drug  delivery  method  in  my  second  year. 

e.  Additionally,  I  have  obtained  the  skills  for  fabricating  the  biodegradable  fibers  which  will  be 
necessary  for  the  second-year  study  of  this  project,  (this  was  originally  scheduled  for  months  17- 
18) 

-  Fundamentals  of  NIRS  system. 

NIR  light  (700  to  900  nm),  where  photon  transport  in  tissue  is  dominated  by  light  scattering  rather  than 
absorption,  has  the  maximal  penetration  depth  (several  centimeters)  in  tissue,  therefore  enabling  in  vivo 
sampling  of  large  tissue  volumes,  such  organs  as  the  breast,  the  brain,  and  skeletal  muscles,  or  deep 
tumors.  Our  multi-channel  NIR  spectroscopy  system  has  the  ability  to  measure  light  signals  from  eight 
different  locations,  but  due  to  the  finite  tumor  size,  we  used  only  four  detectors  to  monitor  tumor  vascular 
oxygenation  dynamics  during  respiratory  challenges.  In  common  with  our  previous  work,9 12  we  assume 
that  oxyhemoglobin  (Hb02)  and  deoxyhemoglobin  (Hb)  are  the  only  significant  absorbing  materials  in 
the  blood-perfused  tumor  tissue.  The  absorption  coefficients  comprise  the  extinction  coefficients  for 
deoxyhemoglobin  and  oxyhemoglobin  multiplied  by  their  respective  concentrations  (Eqs.  1  and  2)  at  730 
nm  and  850  nm,  that  have  been  employed  in  our  multi-channel  NIRS  system. 

Ha™  =  EHb730[Hb]  +  £HbO2730[HbO2],  (1) 

!4850  =  eHb850[Hb]  +  £HbO2850[HbO2] .  (2) 

Based  on  modified  Beer-Lamberf  s  law,  the  data  presented  in  this  report  were  analyzed  using  amplitude 

values  to  find  the  changes  in  light  absorption  (Eq.  3).  By  manipulating  Equations  1-3,  changes  in 

oxygenated  hemoglobin,  deoxygenated  hemoglobin  and  total  hemoglobin  (Hblotai)  concentrations  were 
calculated  from  the  transmitted  amplitude  of  light  through  the  tumor  (Eqs.  4,  5  and  6). 


HaB  -  MaT  =  log  (Ag/AT)  /  L,  (3) 

A[Hb02]  =  [-0.674*log  (AB/\)m  +1.1 17*log  (AB/AT)85°]  /  L,  (4) 

A[Hb]  =  [0.994* log  -  0.376*log  (AB/AT)85°)  /  L,  (5) 

A[Hb](olal  =  A[Hb]  +  A[HbO2]=[0.32*log  (AB/AT)73°  +  0.741  *log  (A^A,)850]  /  L,  (6) 


where  AB  =  baseline  amplitude;  AT  =  transition  amplitude;  L  =  optical  pathlength  between  the  source  and 
detector.  The  constants  contained  in  these  equations  were  computed  with  the  extinction  coefficients  for 
oxy-  and  deoxy-hemoglobin  at  the  two  wavelengths  used.11  In  principle,  L  should  be  equal  to  the  source- 
detector  separation,  d,  multiplied  by  a  differential  pathlength  factor  (DPF),  i.e.,  L=d*DPF.  However,  little 
is  known  about  DPF  for  tumors,  although  a  DPF  value  of  2.5  has  been  used  by  others."  Since  our  focus 
is  on  dynamic  changes  and  relative  values  of  tumor  [Hb02]  in  response  to  oxygen  intervention,  we  have 
taken  the  approach  of  including  the  DPF  in  the  unit,  and  eq.  (4)  becomes  as  follows: 

A[Hb02]  =  [-0.674*log  (AB/AT)m  +1.1 17*log  (Ag/A,.)850]  /  d,  (7) 

where  d  is  the  direct  source-detector  separation  in  cm,  and  the  unit  of  A[Hb02]  in  Eq.  (7)  is  mM/DPF. 
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Accomplishments  during  the  second  period  (months  8-12): 

a.  Dynamic  responses  of  breast  tumor  A[Hb02I  to  oxygen  intervention  monitored  by  multi-channel 
N1RS 

The  dynamic  response  of  tumor  oxygenation  to  oxygen  gas  intervention  has  been  observed  by  multi¬ 
channel  NIRS  system.  In  this  project,  four  detectors  were  placed  on  the  side  of  tumor  and  light  source 
probe  was  located  centrally  on  the  top  of  tumor,  providing  us  with  four  different  regional  signals  of  tumor 
vascular  oxygenation  during  oxygen  inhalation  (Figure  1).  Changes  of  [Hb02],  [Hb],  and  [Hbtoiai]  were 
simultaneously  measured  on  a  group  of  5  rats  during  a  sequence  of  air  to  oxygen  and  then  to  air  for 
studying  the  dynamic  response  of  tumor  oxygenation  to  oxygen  intervention. 


Light  Source 


Air  II  Oxygen 


Isoflu- 
rane 


Multi¬ 
channel 
NIRS  System 


Probe 


Pulse  Oximeter 


Inhaled  Gas 


(b) 


Top  View 


Figure  1.  Schematic  diagram  of  experimental  setup  for  the 
multi-channel  NIR  spectroscopy  system,  (a)  temperature  probe 
was  located  rectally,  and  pulse  oximeter  was  placed  on  the  hind 
paw  to  monitor  the  arterial  oxygen  saturation,  (b)  light  source 
probe  was  placed  on  the  top  center  of  the  tumor,  and  four 
detectors  were  attached  on  four  sides  of  the  tumor  to  detect  the 
regional  signals  during  oxygen  intervention. 


Figure  2.  Time  course  of  tumor  vascular  A[Hb02],  A[Hb], 
and  A[Hbtotai]  for  a  representative  13762NF  breast  tumor 
with  the  inhaled  gas  under  the  sequence  of  air-oxygen-air. 
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A  representative  result  is  shown  in  Fig.  2.  After  10  minutes  of  air  breathing  measurement  as  the 
baseline,  the  inhaled  gas  was  switched  from  air  to  oxygen,  causing  a  sharp  increase  in  A[Hb02]  (p  < 
0.0001  1  min.  after  gas  switch)  followed  by  a  further  gradual,  but  significant  increase  over  the  next  13 
min  (P  <  0.0001).  However,  A[Hbtotai]  first  dropped  rapidly  after  switching  gas  from  air  to  oxygen  and 
then  became  stabilized  during  the  next  13  minutes  of  oxygen  inhalation.  This  behavior  is  expected  since 
oxygen  intervention  may  cause  vasoconstriction  of  artery  and  result  in  blood  flow  decreases,  which 
reduces  blood  volume.  After  returning  to  air  inhalation,  A[Hb02]  decreased  and  A[Hbtotai]  increased,  both 
reaching  approximately  their  respective  baseline  level.  This  result  clearly  shows  that  an  NIRS  can 
monitor  the  changes  of  blood  oxygenation  or  blood  volume  in  tumor  dynamically  and  in  real-time. 

b.  Body  weight  and  tumor  volume  changes  during  chemotherapy 

The  tumor  type,  which  we  are  using,  is  murine  mammary  adenocarcinomas  13762NF,  and  this  tumor  type 
is  highly  responsive  to  alkylating  agents  and  platinum  chemotherapeutic  agents.14  Cyclophosphamide, 
one  of  alkylating  agents,  has  been  used  in  various  types  of  cancer  such  as  leukemia,  ovarian  cancer  and 
breast  cancer.15  We  have  injected  cyclophosphamide  into  the  rat  through  i.  p.  with  a  single  high  dose  of 
200  mg/kg  for  this  experiment  to  see  whether  we  can  detect  acute  effects  of  cyclophosphamide  on  tumor 
vascular  oxygenation  dynamics.  Some  rats  were  failed  to  survive  during  the  treatment,  but  the  others 
survived  with  apparent  reduction  of  tumor  volume  during  the  10  days  after  i.  p.  injection  of 
cyclophosphamide.  Figures  3a  and  3b  show  changes  in  rat  body  weight  and  tumor  volume,  respectively, 
before  and  after  administration  of  cyclophosphamide.  (Day  0  was  the  day  for  the  injection.)  These  Figures 
show  clearly  decreases  in  rat  body  weight  as  well  as  in  tumor  volume  for  5  experimental  rats. 


Day  Day 


(a)  (b) 

Figure  3.  Changes  of  rat  body  weight  (a)  and  tumor  volume  (b)  during  chemotherapy. 

c.  Intratumoral  heterogeneity  of  vascular  oxygenation  observed  by  multi-channel  NIRS 

Solid  tumors  are  known  to  exhibit  heterogeneous  blood  flow  distribution. 16,1  Intensive  studies  from 
Mason’s  group  using  l9F  MR  [FREDOM  (Fluorocarbon  Relaxometry  using  Echo  planar  imaging  for 
Dynamic  Oxygen  Mapping)]  have  revealed  intratumoral  heterogeneity  of  p02  distribution  and  also 
heterogeneous  response  to  hyperoxic  gas  breathing.'819,20  The  findings  of  severe  p02  heterogeneity  in 
tumors  suggest  the  heterogeneous  distribution  of  blood  flow,  since  tissue  p02  level  is  determined  by  a 
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balance  between  the  supply  of  oxygen  from  blood  vessels  and  the  oxygen  consumption  rate  of  tissue 
cells.21  Unlike  FREDOM,  near-infrared  spectroscopy  (NIRS)  techniques  measure  hemoglobin  oxygen 
saturation  and  concentration  in  vivo,  providing  information  on  vascular  oxygenation  of  the  measured 
sample/organ,  such  as  the  brain3,4'23  and  animal  tumors. 9,1123  Therefore,  multi-channel  NIRS  can  detect 
signals  from  many  locations  of  a  solid  tumor,  enabling  us  to  study  vascular  heterogeneity  of  a  tumor. 

In  our  experiment,  we  have  applied  four  detectors  and  one  light  source  on  the  rat  breast  tumors,  and  the 
setup  has  been  shown  in  Fig.  1.  After  10  minutes  of  baseline  measurement  with  air  breathing,  gas  was 
switched  to  pure  oxygen,  which  caused  a  rapid  increase  of  tumor  A[Hb02].  These  changes  were  measured 
simultaneously  from  4  locations  of  the  tumor,  and  Fig.  4  shows  a  representative  set  of  data.  Red  circles 
show  the  raw  data  measured  by  multi-channel  NIRS,  and  the  solid  black  lines  represent  the  fitted  curves 
using  our  double-exponential  model  for  hemodynamics  during  oxygen  intervention9.  As  we  can  see  from 
these  figures,  none  of  the  data  from  four  detectors  are  identical  even  though  some  are  showing  quite 
similar  trends  of  A[Hb02]. 
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Figure  4.  Dynamic  changes  of  tumor  [HbOi]  from  four  detectors  in  a  rat  breast  tumor  (Tumor  A  with  a  volume  of  8.2  cm  ). 
The  rising  part  of  A[Hb02]  from  the  four  detectors  was  fitted  using  double-exponential  expressions  and  Fig.  4(a)-4(d)  are  taken 
from  detectors  #l-#4,  respectively.  In  this  case,  the  tumor  was  not  treated  yet. 
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Table  1.  Summary  of  vascular  oxygen  dynamics  determined  at  the  four  detectors  from  Figure  4  for  Tumor  A. 


Detector 

A1 

A2 

xl 

t2 

A1/A2  (=Yi/Y2) 

t1/t2 

/.//2 

No.  1 

0.013 

0.027 

0.48 

24.3 

0.48 

0.02 

24 

No.  2 

0.029 

0.026 

0.09 

29.8 

1.12 

0.003 

373 

No.  3 

0.008 

0.036 

0.15 

10.6 

0.22 

0.014 

16 

No.  4 

0.021 

0.037 

0.12 

7.25 

0.57 

0.017 

34 

To  compare  the  data  taken  from  four  locations  of  the  tumor  more  clearly,  the  time  constants  and 
amplitudes  from  the  four  fitted  curves  are  summarized  in  Table  1.  The  ratios  of  Y1/Y2  and  /j//2  basically 
characterize  tumor  vascular  structure  and  blood  perfusion  within  the  volume  of  tumor  interrogated  by 
light.9  In  principle,  when  Y1/Y2  is  close  to  1,  it  implies  that  the  measured  optical  signal  results  equally  from 
both  region  1  (i.e.,  well  perfused  region)  and  2  (i.e.,  poorly  perfused  region);  if  Y1/Y2  <  1,  the  measured 
signal  results  more  from  region  2  than  region  1  [Figures  3(a),  3(c)  and  3(d)].  As  Table  1  demonstrates, 
only  detector  #2  has  a  ratio  of  Y1/Y2  slightly  higher  than  1,  and  the  readings  from  detectors  #1,  #3  and  #4 
have  the  ratios  of  Y1/Y2  less  than  1.  This  may  suggest  that  the  tumor  volume  that  was  optically  interrogated 
by  detector  #2  consists  of  more  volumes  from  well  perfused  region  than  poorly  perfused  region,  while 
most  of  other  tumor  volumes  detected  by  detectors  #1,  #3  and  #4  are  composed  of  more  poorly  perfused 
regions.  Furthermore,  all  the  ratios  of  f\lf2  from  Tumor  A  at  the  four  locations  are  much  greater  than  1, 
indicating  that  the  blood  perfusion  rate  in  well  perfused  region  is  much  greater  than  that  in  poorly 
perfused  region.  Especially,  f\lf2  from  detector  #2  is  10  to  20  times  higher  than  those  from  detectors  #1, 
#3,  and  #4,  showing  a  high  level  of  intratumoral  heterogeneity  in  dynamic  vascular  structure. 


d.  Monitoring  vascular  hemodynamics  of  breast  tumors  before  and  after  chemotherapy 

The  hemodynamics  of  tumor  during  oxygen  intervention  was  measured  before  and  after  administration  of 
cyclophosphamide.  The  representative  data  from  tumor  B  are  shown  in  Fig.  5.  As  it  was  mentioned  before, 
there  were  4  detectors  on  the  surface  of  tumor,  and  Fig.  5  shows  the  acute  and  then  gradual  changes  of 
[Hb02]  after  switching  the  breathing  gas  from  air  to  oxygen,  and  the  data  were  observed  by  detector  #1 
from  day  0  to  day  5.  This  figure  clearly  demonstrates  that  we  can  observe  significant  changes  in  tumor 
hemodynamics  after  chemotherapy  by  giving  respiratory  challenges.  The  fitted  parameters  from  our 
mathematical  model  are  summarized  in  Table  2  to  compare  the  changes  in  hemodynamic  parameters 
before  and  after  administration  of  cyclophosphamide.  The  rising  part  of  A[Hb02]  from  detector  #1  was 
fitted  using  either  single-  or  double-  exponential  expression. 


Table  2.  Summary  of  vascular  oxygen  dynamics  determined  at  detector  #1  from  Tumor  B  before  and  after 

cyclophosphamide  administration. 


Day 

A1 

A2 

xl 

t2 

A1/A2(^y,/y2) 

t1/t2 

Mi 

0 

0.044 

0.031 

0.23 

25.21 

1.42 

0.0091 

156 

1 

0.032 

0.033 

0.13 

11.36 

0.97 

0.0114 

85 

3 

0.0087 

0.014 

0.089 

8.36 

0.62 

0.0106 

58 

5 

0.0087 

1.27 

10 


Time(min) 

Figure  5.  Dynamic  changes  of  [Hb02]  taken  at  detector  #1  from  a  rat  breast  tumor  before  and  after  administration  of 
cyclophosphamide  (Tumor  B  with  a  volume  of  5.3  cm3  at  day  0).  The  rising  part  of  A[Hb02]  from  detector  #1  was  fitted  using 
either  single-  or  double-exponential  expression. 


At  day  0,  we  can  see  that  Yt/y2  is  higher  than  1,  indicating  that  the  measured  signal  results  more  from  the 
well  perfused  region  than  poorly  perfused  region.  However,  this  ratio  becomes  less  than  1  after  injection 
of  cyclophosphamide  (Day  1  and  3).  This  may  be  explained  by  destruction  of  vascular  structure  in  tumor 
after  chemotherapy  as  follows.  After  administration  of  cyclophosphamide,  it  will  be  circulated  in  blood 
vessels  and  will  be  more  delivered  to  the  tumor  cells  in  the  well  perfused  region  than  those  in  the  poorly 
perfused  region.  This  will  cause  more  apoptosis  of  tumor  cells  in  the  well  perfused  region  than  that  in  the 
poorly  perfused  region,  eventually  causing  decreases  in  tumor  volume  in  the  well  perfused  region  more 
than  that  in  the  poorly  perfused  region.  Then,  the  tumor  volume  optically  interrogated  by  detector  #1  will 
have  less  contribution  from  the  well  perfused  region  than  the  poorly  perfused  region,  leading  to  a 
decrease  of  Yi/y2,  after  administration  of  cyclophosphamide. 

After  cyclophosphamide  administration,  the  perfusion  rate  ratio,  f}/f2,  was  also  decreased.  At  day 
is  very  high,  which  means  there  is  a  big  difference  of  perfusion  rate  between  the  well  perfused  and  poorly 
perfused  region  in  tumor.  However,  this  ratio  significantly  decreased  at  day  1  and  3  after 
cyclophosphamide  administration,  representing  that  the  perfusion  rate  gap  between  the  well  perfused 
region  and  poorly  perfused  region  became  much  smaller  than  that  from  day  0.  At  day  5,  the  changes  of 
[Hb02]  during  oxygen  intervention  does  not  show  any  bi-phasic  behavior  anymore,  and  it  was  fitted  by  a 
single-exponential  model.  This  may  indicate  that  most  of  tumor  cells  in  the  well  perfused  region  possibly 
destroyed  by  the  effect  of  cyclophosphamide,  and  vascular  structure  has  been  totally  changed. 

Similar  experiments  have  been  performed  on  more  rat  tumors  using  the  same  methodology  and  analysis. 
A  decrease  of  maximum  A[Hb02]  (A[Hb02]max)  during  oxygen  inhalation  was  clearly  observed  from  all 
tumors  after  cyclophosphamide  administration.  Figure  6  shows  a  set  of  dynamic  changes  of  [Hb02]  right 
after  switching  the  gas  from  air  to  oxygen  at  day  0  and  3  after  cyclophosphamide  administration.  All  of 
the  four  tumors  showed  rapid  increase  of  A[Hb02]  to  the  oxygen  inhalation,  but  A[Hb02]max  was  reduced 
at  day  3  by  20  to  70%  compared  to  A[Hb02]„m  at  day  0.  We  notice  that  the  decrease  of  A[Hb02]nBX  at  day 
3  is  mainly  from  the  decrease  of  A1  (amplitude  of  rapid  increase  in  A[Hb02])  rather  than  A2  (amplitude 
of  gradual  increase  in  A[Hb02]).  Table  3  summarizes  the  fitted  parameters  of  Fig.  6.  As  seen  from  Table 
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3,  Yi/y2*s  decreased  for  all  four  tumors  at  day  3  compared  to  day  0.  The  perfusion  rate  ratio,  f\lf2,  has  also 
been  decreased  for  all  tumors  at  day  3  compared  to  day  0.  These  observations  are  in  good  agreement  with 
the  result  shown  in  Figure  5  and  Table  2.  All  of  the  results  may  suggest  that  the  administration  of 
cyclophosphamide  treats  the  tumor  well  with  well  perfused  vasculature,  but  not  well  for  the  tumors  with 
poorly  perfused  vasculature.  Furthermore,  our  study  shown  here  proves  that  we  can  monitor  the  effects  of 
chemotherapy  in  tumor  by  detecting  changes  of  hemodynamics  during  respiratory  challenges  using  near 
infrared  spectroscopy. 


Figure  6.  Dynamic  changes  of  [Hb02]  after  switching  the  gas  from  air  to  oxygen  from  four  individual  rat  breast  tumors.  Red 
open  circles  and  solid  lines  indicate  raw  data  of  the  changes  of  [Hb02]  at  day  0  and  fitted  curve,  respectively.  Blue  open  squares 
and  solid  lines  indicate  raw  data  of  the  changes  of  [Hb02]  at  day  3  after  cyclophosphamide  administration  and  fitted  curve, 

respectively. 
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Table  3.  Summary  of  vascular  oxygen  dynamics  determined  at  day  0  and  day  3  from  four  tumors. 


Tumor 

Day 

A1 

A2 

xl 

t2 

A1/A2 

(=Yl/Y2) 

xl/  t2 

Mi 

Fig.6a 

0 

0.034 

0.015 

0.34 

23.92 

2.27 

0.014 

159 

3 

0.013 

0.016 

0.11 

3.57 

0.81 

0.031 

26 

Fig.6b 

0 

0.021 

0.036 

0.24 

20.51 

0.58 

0.012 

50 

3 

0.005 

0.045 

0.14 

34.10 

0.11 

0.004 

27 

Fig.6c 

0 

0.040 

0.046 

0.17 

17.87 

0.87 

0.010 

91 

3 

0.015 

0.081 

0.07 

29.54 

0.19 

0.002 

76 

Fig.6d 

0 

0.043 

0.066 

0.25 

49.33 

0.65 

0.005 

129 

3 

0.0053 

0.019 

0.10 

12.23 

0.28 

0.008 

34 

Key  Research  Accomplishments: 

a.  Obtained  the  main  knowledge  and  skills  for  conducting  this  project,  (training  in  handling  laboratory 
animals,  studying  the  fundamentals  and  operation  in  the  NIR  oximeter  and  19F  MR  p02  mapping 
method,  studying  the  drug  delivery  and  also  tumor  biology,  learning  the  fabrication  of  biodegradable 
fiber.) 

b.  Observed  the  changes  of  tumor  volume  and  body  weight  of  rats  during  chemotherapy.  This  clearly 
showed  the  response  of  tumor  to  the  chemotherapy. 

c.  Studied  the  intratumoral  heterogeneity  of  breast  tumors  by  multi-channel  NIRS.  Changes  of  [Hb02], 
[Hb]  and  [Hb,ota|]  during  gas  interventions  were  monitored  by  multi-channel  NIRS  in  real  time,  and 
the  data  clearly  revealed  that  there  is  dynamic  intratumoral  heterogeneity  in  tumor  vascular  structure 
during  oxygen  intervention. 

d.  Investigated  the  dynamic  response  of  tumor  vascular  A[HbC>2]  to  oxygen  intervention  before  and 
after  cyclophosphamide  administration.  The  ratios  of  Y1/Y2  and  /1//2  have  been  changed  during 
therapy,  indicating  that  changes  in  hemodynamics  of  AfHbCL]  during  oxygen  intervention  before 
and  after  cyclophosphamide  can  be  applied  to  find  the  efficacy  of  chemotherapeutic  agent. 

Reportable  Outcomes 

Manuscripts  for  Peer-reviewed  Journals: 

1)  Yueqing  Gu,  Vincent  A.  Bourke,  Jae  G.  Kim,  Anca  Constantinescu,  Ralph  P.  Mason,  Hanli  Liu, 
“Dynamic  response  of  breast  tumor  oxygenation  to  hyperoxic  respiratory  challenge  monitored  with  three 
oxygen-sensitive  parameters”,  Applied  Optics,  42(16),  2960-2967  (2003) 

2)  Hanli  Liu,  Yueqing  Gu,  Jae  G.  Kim,  Ralph  P.  Mason,  “Near  infrared  spectroscopy  and  imaging  of 
tumor  vascular  oxygenation”.  Imaging  in  Biological  Research,  Part  B  (P.  Machael  Conn  ed.)  Methods  in 
Enzymology,  386,  pp.  349-378  (2004) 

3)  Jae  G.  Kim,  Hanli  Liu,  “Investigation  of  Breast  Tumor  Hemodynamics:  Dynamic  Phantoms  and  FEM 
Simulations”,  in  preparation  and  will  be  submitted  to  applied  optics. 

4)  Jae  G.  Kim,  Anca  Constantinescu,  Ralph  P.  Mason,  Hanli  Liu,  “Non-Uniform  Tumor  Vascular 
Oxygen  Dynamics  Monitored  by  Multi-Channel  Near-Infrared  Spectroscopy,  in  preparation  and  will  be 
submitted  to  J.  Biomedical  Optics. 
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Presentations  and  Proceeding  papers: 

l)  Jae  G.  Kim,  and  Hanli  Liu,  “Investigation  of  breast  tumor  hemodynamics  using  tumor  vascular 
phantoms  and  FEM  simulations”,  in  Biomedical  Topical  Meetings  on  CD-ROM  (The  Optical  Society  of 
America,  Washington,  DC,  2004),  WF16. 


Conclusions 

From  the  work  that  I  have  conducted  up-to-date,  the  following  conclusions  can  be  drawn. 

a.  The  dynamic  responses  of  A[Hb02]  after  switching  the  breathing  gas  from  air  to  oxygen  proved 
that  there  exists  dynamic  intratumoral  heterogeneity  in  breast  tumor  vascular  structure  during 
oxygen  intervention.  During  the  oxygen  inhalation,  all  four  detectors  showed  rapid  increases  of 
A[Hb02]  followed  by  a  gradual  increases.  (The  bi-phasic  behavior  of  A[Hb02]  implies  that  the 
signal  in  the  fast  phase  results  from  well  perfused  regions,  and  the  signal  in  the  slow  phase  results 
from  poorly  perfused  regions.)  However,  two  ratios  (Y1/Y2  and  /t//2)  from  the  double-exponential 
fitting  to  the  data  showed  large  differences  among  the  signals  taken  from  each  detector.  The 
differences  suggest  that  some  parts  of  tumor  consist  of  more  well  perfused  regions  than  other 
parts.  The  results  demonstrate  the  capability  of  NIR  spectroscopy  as  a  noninvasive,  real  time 
monitoring  tool  for  tumor  vascular  oxygenation. 

b.  Changes  in  hemodynamic  parameters  before  and  after  chemotherapy  were  monitored  from  breast 
rat  tumors.  The  results  showed  that  A[Hb02]maX  during  oxygen  inhalation  has  been  decreased 
largely  after  cyclophosphamide  administration,  which  is  consistent  with  the  reduction  of  tumor 
size.  Comparing  the  two  ratios  (Y1/Y2  and  /i//2)  before  and  after  the  drug  injection  revealed  that 
tumor  blood  volume  in  the  well  perfused  region  may  be  decreased  faster  than  that  in  the  poorly 
perfused  region.  Our  results  may  imply  that  destruction  of  tumor  cells  occurred  more  efficiently 
in  the  well  perfused  region  than  in  the  poorly  perfused  region  of  tumor.  This  is  expected  because 
the  administered  drug  will  be  delivered  much  easier  to  the  tumor  cells  near  the  well  perfused 
region  than  those  far  from  the  blood  vessels.  All  of  our  data  obtained  in  Year  1  prove  that 
measuring  hemodynamics  of  breast  tumor  during  gas  intervention  can  be  feasibly  used  for 
monitoring  the  efficacy  of  cancer  therapy,  such  as  radiation  therapy,  photodynamic  therapy,  and 
chemotherapy. 
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ABSTRACT 


In  a  previous  study,  switching  gas  from  air  to  carbogen  caused  a  sharp  increase  of  oxyhemoglobin 
concentration  in  breast  tumors  followed  by  a  further  gradual,  but  significant,  increase.  This  bi-phasic 
behavior  of  tumor  vascular  oxygenation  was  fitted  using  a  mathematical  model  based  on  the  hypothesis 
that  the  measured  signal  results  from  both  well-perfused  and  poorly  perfused  region  in  tumor.  To  prove 
our  hypothesis,  we  have  designed  and  fabricated  a  breast  tumor  dynamic  phantom  in  order  to  perform 
three  channel  experiments  with  the  multi-channel  NIR  system.  Tumor  dynamic  phantoms  are  made  of 
gel  that  encapsulated  transparent  tubing  through  which  dye  flow  rates  were  controlled  via  two  syringe 
pumps.  Tumor  phantom  experiments  showed  that  two  time  constants  in  tumor  hemodynamics  are  indeed 
from  two  different  perfusion  rates.  We  also  conducted  computer  simulations  to  further  support  our  tumor 
phantom  measurements  and  to  understand  the  effects  of  tumor  heterogeneity  on  the  multi-channel  NIR 
readings  using  a  finite  element  method.  (FEM)  The  FEM  simulation  results  were  well  matched  with  the 
results  from  tumor  phantom  experiments  and  thus  these  two  methods  proved  our  hypothesis.  This  study 
demonstrates  that  by  having  an  appropriate  model,  we  can  monitor  the  therapeutic  effects  of  cancer 
treatments  using  non-invasive  NIRS  and  it  may  have  prognostic  value,  providing  insight  into  tumor 
vascular  development  and  angiogenesis. 

Key  Terms:  Tumor  vascular  oxygenation,  Finite  element  method,  NIR  Spectroscopy,  Tumor  phantom, 
Carbogen 
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INTRODUCTION 


Solid  tumors  develop  regions  of  hypoxia  during  their  growth  due  to  an  imbalance  between  the 
rate  of  tumor  cell  proliferation  and  branching  of  the  blood  vessels.1,2,3  Tumor  hypoxia  can  contribute  to 
the  failure  of  radiotherapy,4  some  forms  of  chemotherapy,5  and  photodynamic  therapy.6  In  addition,  a 
number  of  clinical  trials  have  found  that  patient  survival  measured  either  as  tumor  regression  or  as  local 
control  depends  highly  on  tumor  oxygenation.7  Therefore,  increasing  tumor  oxygenation  could  improve 
cancer  therapy.  As  one  means  to  improve  tumor  oxygenation,  breathing  hyperoxic  gas  has  been  used  to 
enhance  treatment  efficacy.  However,  poor  tumor  perfusion  may  prohibit  the  oxygen  delivery  from  blood 
vessels  to  tumor  tissue  cells.  Therefore,  measurements  of  local  vascular  oxygenation  and  distribution  of 
tumor  perfusion  could  be  important  for  tumor  treatment  planning  and  to  assess  methods  designed  to 
modulate  tumor  oxygenation. 

Solid  tumors  are  also  known  to  exhibit  heterogeneous  blood  flow  distribution.8  9  There  are  various 
methods  to  study  tumor  perfusion  heterogeneity,  such  as  Doppler  ultrasound,10  dynamic  contrast  MRI,11 
and  the  use  of  tumors  grown  in  window  chambers.12  Intensive  studies  from  Mason’s  group  using  19F  MR 
[FREDOM  (Fluorocarbon  Relaxometry  using  Echo  planar  imaging  for  Dynamic  Oxygen  Mapping)]  have 
revealed  intratumoral  heterogeneity  of  p02  distribution  and  also  heterogeneous  response  to  hyperoxic  gas 
breathing.  13  ’ 14  ’ 15  The  findings  of  severe  p02  heterogeneity  in  tumors  suggest  the  heterogeneous 
distribution  of  blood  flow,  since  tissue  p02  level  is  determined  by  a  balance  between  the  supply  of  oxygen 
from  blood  vessels  and  the  oxygen  consumption  rate  of  tissue  cells.16  Unlike  FREDOM ,  near-infrared 
spectroscopy  (NIRS)  techniques  measure  hemoglobin  saturation  and  concentration  in  vivo ,  providing 
information  on  vascular  oxygenation  of  the  measured  sample/organ.  NIRS  has  been  widely  applied  to 
investigate  hemoglobin  oxygenations  of  muscles,  *  ’  the  brain,"  ’  ’  and  animal  tumors.  *  ’ 

Our  previous  in  vivo  animal  study  has  clearly  demonstrated  that  carbogen  (95%  C02  and  5%  02) 
inhalation  can  improve  the  vascular  oxygen  level  of  breast  tumor.  The  observed  changes  in  tumor 
vascular  oxygenation  exhibited  a  rapid  increase,  followed  by  a  gradual  but  significant  increase,  in 
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response  to  carbogen  intervention.  (Figure  1)  To  explain  these  biphasic  behavior  of  tumor  hemodynamics, 
we  established  a  mathematical  model  based  on  the  Kety’s  approach26  and  we  formed  a  hypothesis  that 
tumor  vasculature  is  comprised  with  a  well-perfused  and  poorly  perfused  region  to  explain  why  there  are 
two  different  time  constants  in  the  A[Hb02]  data.  The  model  further  allows  us  to  associate  the  signal 
amplitudes  to  the  ratio  of  vascular  coefficients  and  the  ratio  of  the  perfusion  rates  in  the  two  different 
regions. 

Even  though  we  associated  our  fitting  parameters  with  physiological  factors  such  as  vascular 
coefficients  and  perfusion  rates  mathematically,  it  was  not  proved  at  that  time.  To  provide  actual  support 
for  this  model  and  to  further  investigate  heterogeneity  of  breast  tumor  vasculature,  we  designed  and 
developed  a  breast  tumor  hemodynamic  phantom  in  order  to  perform  three  channel  experiments  with  the 
multi-channel  NIRS  system.  To  further  support  our  tumor  hemodynamic  phantom  study,  we  used 
FEMLAB  software  to  simulate  the  dynamic  process  of  light  distribution  as  well  as  the  NIR  measurement 
in  breast  tumors. 

MATERIALS  AND  METHODS 

Mathematical  Model  of  Tumor  Vascular  Oxygenation 

In  our  previous  report,25  we  followed  an  approach  used  to  measure  regional  cerebral  blood  flow 
(rCBF)  with  diffusible  radiotracers,  as  originally  developed  by  Kety26  in  the  1950’s.  By  applying  Fick’s 
principle  and  defining  y as  the  ratio  of  Hb02  concentration  change  in  the  vascular  bed  to  that  in  veins,  we 
arrived  at  Eq.  (1): 

AHbO2vascu,ature(0  =  yHJ  l-exp(-yr/y)]  =  A[  l-exp(-r/x)]  (1) 
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where  y  is  the  vasculature  coefficient  of  the  tumor  (=A[Hb02]vasculalure/A[Hb02]vem),  Ho  is  the 
arterial  oxygenation  input  and  /is  the  blood  perfusion  rate. 

If  a  tumor  has  two  distinct  perfusion  regions  and  the  measured  signal  results  from  both  regions 
(Figure  2),  then  it  is  reasonable  to  include  two  different  blood  perfusion  rates,/  and/,  and  two  different 
vasculature  coefficients,  y,  and  y>,  in  the  model.  Therefore,  Eq.  (1)  can  be  modified  to  count  for  the 
double  exponential  feature  observed  in  the  experiments: 

AHbO2vasculaturc(0  =  KH„[l-exp(-/r/n)]  +  ^HJl-exp/Zt/ft)] 

=  Ai[l-exp(-//Ti)]  +  A2[l-exp(-r/t2)]  (2) 

where/i  and  Y\  are  the  blood  perfusion  rate  and  the  vasculature  coefficient  in  region  1  for  the  well 
perfused  region,  respectively;  /2  and  Yi  have  the  same  respective  meanings  in  region  2  for  the  poorly 
perfused  region,  and  Ai  =  ftHo,  A2  =  ^Ho,  X\  =  Y\/fu  t2  =  Yitfi-  Then,  if  Ai,  A2,  Ti,  and  t2  are  determined 
by  fitting  the  measurements  with  the  model,  we  can  obtain  the  ratios  of  two  vasculature  coefficients  and 
the  two  blood  perfusion  rates: 

Ki  _  A  f\  _  AM 2  (3) 

Yi  f2 

With  these  two  ratios,  we  are  able  to  understand  more  about  tumor  physiology,  such  as  tumor 
vasculature  structure  and  blood  perfusion.  In  this  report,  we  used  three-channel  NIR  spectroscopy  to 
evaluate  our  tumor  hemodynamics  model  by  quantifying  A,  x,  Y1/Y2  and/i//2  at  three  different  locations  on 
the  tumor  phantoms. 

Tumor  Hemodynamic  Phantom  Design  and  Measurements 
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To  represent  vascular  blood  vessels  in  tumor,  we  designed  a  vascular  modeling  device  (VMD)  by 
winding  a  small  diameter  tube  around  a  big  diameter  core  tube.  An  example  of  such  a  VMD  is  shown  in 
Figure  3.  VMD-1  was  fabricated  by  wrapping  ethyl  vinyl  acetate  microbore  tubing  (0.51  mm  ID)  around 
tygon  lab  tubing  (14.4  mm  OD),  and  VMD-2  was  fabricated  by  winding  tygon  silicone  tubing  (0.80  mm 
ID)  around  tygon  lab  tubing  (14.4  mm  OD)  to  represent  small  diameter  blood  vessel  and  large  diameter 
blood  vessel,  respectively.  All  tubing  materials  were  purchased  from  Cole-Parmer  Company  (Vemon 
Hills,  IL). 

Tumor  dynamic  phantom  was  fabricated  by  embedding  two  VMDs  into  a  cylinder  shape  of  gel, 
which  represents  tissues.  Tissue  mimic  gel  was  made  by  mixing  50  g  of  gel  powder  (Sigma,  Gelatin  Type 
A,  St.  Louis,  MO)  into  350  ml  of  boiling  water,  and  solution  was  stirred  thoroughly  until  gel  powder  was 
dissolved  completely.  When  the  solution  was  cooled  down  to  around  50  °C,  200  ml  of  intralipid 
(Intralipid"  20%,  Baxter  Healthcare  Corp.,  Deerfield,  EL)  was  added  and  mixed  thoroughly  to  simulate 
light  scattering  in  tumor  tissues.  Just  before  solution  started  to  be  solidified,  it  was  poured  into  a 
cylindrical  container  (diameter=4.5  cm,  height=3.5  cm),  which  contained  two  VMDs  inside.  After 
completely  being  cooled  down,  the  solution  became  a  soft  gelatin  phantom  with  the  two  VMDs 
embedded  inside.  The  optical  properties  of  tumor  gel  phantom  were  measured  using  near  infrared  tissue 
spectrometer  (model:96208,  ISS  Inc.,  Champaign,  IL)  and  those  values  were  close  to  tissue  optical 
properties  (na=  0.032  cm'1,  ps’=9.2  cm'1  at  750  nm). 

For  our  experiments,  two  tumor  hemodynamic  phantoms  were  fabricated.  Phantom  1  contained 
one  VMD-1  and  one  VMD-2  while  two  VMD- Is  were  embedded  into  a  cylindrical  gel  for  phantom  2.  A 
stream  of  dye  solution  (diluted  black  ink,  pa=1.5  cm'1  at  730  nm)  was  used  to  go  through  the  VMDs  to 
simulate  a  blood  flow  through  tumor  vasculature.  By  applying  same  dye  flow  rate  into  VMD-1  and 
VMD-2  in  phantom  1,  we  would  see  the  effects  of  blood  vessel  length  or  diameter  on  bi-phasic  behavior 
of  tumor  blood  oxygenation  during  carbogen  inhalation.  Meanwhile,  injecting  dye  into  two  VMD-ls  in 
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phantom  2  with  different  flow  rates  would  enable  us  to  find  how  blood  flow  rate  affects  tumor’s  bi-phasic 
increase  of  oxyhemoglobin  during  carbogen  inhalation.  In  this  way,  we  could  mimic  a  hemodynamic 
process  of  a  breast  tumor  by  changing  the  diameter  or  the  length  of  the  small  tube  or  by  controlling  the 
flow  rate  of  the  dye  solution. 

Schematic  experimental  setup  for  tumor  hemodynamic  phantom  study  is  shown  in  Figure  4a.  As 
shown  here,  NIR  light  (730  nm)  from  multi  channel  NIRS  was  delivered  to  the  tumor  hemodynamic 
phantom,  and  three  optical  detectors  placed  on  the  side  of  cylindrical  phantom  collected  signals,  which 
would  be  processed  in  the  computer.  Dye  was  injected  into  VMDs  using  B-D™  disposable  syringe 
(Cole-Parmer,  Vernon  Hills,  IL)  by  two  syringe  pumps  (model  KDS200,  KdScientific  Inc.,  New  Hope, 
PA),  and  dye  wastes  were  collected  in  the  waste  beaker.  We  used  two  syringe  pumps  so  that  we  could 
control  dye  flow  rate  for  each  VMD.  Figure  4b  shows  more  detail  setup  of  light  source  and  three  optical 
detectors  around  the  tumor  phantom.  Light  source  was  placed  between  two  VMDs,  and  detector  D3  was 
located  across  the  light  source  in  the  transmission  mode  so  that  it  would  detect  signal  passing  through 
both  VMDs.  Other  two  detectors  (D1  and  D2)  were  placed  in  the  reflectance  geometry  with  respect  to 
light  source  so  that  D1  and  D2  would  get  the  NIR  signals  mostly  from  only  one  of  two  VMDs. 

Experimental  procedures  are  as  follows.  For  the  tumor  phantom  1,  VMD-1  and  VMD-2  were 
filled  with  water  first  as  a  baseline,  and  then  dye  was  injected  into  the  VMD-2  first  with  a  flow  rate  of  20 
ml  hr'1  to  simulate  the  light  absorption  effect  of  blood,  followed  by  a  washout  with  water.  For  the  second 
step,  dye  was  injected  into  the  VMD-1  with  the  same  flow  rate  as  that  in  the  first  step  and  also  followed 
by  a  washout  of  water.  The  last  step  was  introducing  dye  into  both  VMD-1  and  VMD-2  with  the  same 
flow  rate  of  20  ml  hr'1.  By  flowing  dye  into  a  tumor  phantom  1  in  this  way,  we  expected  to  see  two 
exponential  increases  of  light  absorption  changes  because  dye  flow  velocities  will  be  different  between 
VMD-1  and  VMD-2  due  to  different  sizes  of  tube  diameter  even  though  the  applied  dye  flow  rate  was 
same  for  both  VMDs. 

For  the  tumor  phantom  2,  we  first  filled  two  VMD- Is  with  water  and  then  dye  was  injected  into 
the  top  VMD-1  with  a  flow  rate  of  20  ml  hr  1  followed  by  a  washout  with  water.  In  the  second  step,  dye 
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was  injected  into  the  bottom  VMD-1  with  a  20  ml  hr'1  flow  rate  and  was  washed  out  with  water  again. 
Thirdly,  we  introduced  dye  into  both  VMD-ls  with  a  same  flow  rate  of  20  ml  hr  1  following  washout 
water.  Finally,  dye  was  injected  into  two  VMD-ls  with  two  different  flow  rates  (5  and  20  ml  hr  1  for  top 
and  bottom  VMD-1,  respectively)  and  was  washed  out  with  water.  Changes  of  light  absorption  were 
measured  through  the  whole  experiment,  and  the  time  constants  for  the  dynamic  changes  in  light 
absorption  were  obtained  by  fitting  the  data  with  our  mathematical  models  using  Kaleidagraph  software 
(Synergy  Software,  Reading,  PA). 


FEM  Simulation 

To  further  support  our  tumor  hemodynamic  phantom  study,  we  used  FEMLAB  software 
(COMSOL  Inc.  Burlington,  MA),  which  is  based  on  the  finite  element  method,  to  simulate  our  tumor 
phantom  2  experiment.  It  uses  numerical  approaches  to  solve  partial  differential  equations  (PDE)  in 
modeling  and  simulating  various  engineering  problems.  The  geometry  of  FEM  simulation  for  tumor 
hemodynamic  phantom  is  shown  in  the  figure  5  and  it  represents  a  2  dimensional  horizontal  cross  section 
of  our  phantom  2.  El  represents  the  body  of  tumor  phantom  (diameter  =  4  cm),  and  R5  shows  the 
location  of  the  light  source.  Furthermore,  each  tube  used  in  the  dynamic  phantom  is  represented  by  two 
rectangles  (0.1  x  2.4  cm  rectangle)  as  a  blood  vessel  network:  R1  and  R2  are  representing  vessels  with  a 
fast  flow  rate,  and  R3  and  R4  denote  vessels  with  a  slow  flow  rate. 

Equation  4  is  the  diffusion  equation  that  has  been  applied  to  predict  photon  distribution  in 
phantom  as  dye  flows  through  two  VMDs. 

(l / c\d / DV2</>(r,t)  +  jun<p(r,t)  =  S(r,t)  (4) 
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where  <p( r,t)  is  the  diffuse  photon  fluence  rate  at  the  position  r,  c  is  the  speed  of  light  in  the  tissue,  S(r,t) 
describes  the  photon  source,  and  D  =  [3(pa  +  p/)]*1  is  the  photon  diffusion  coefficient. 

For  the  boundary  conditions,  the  partial  current  boundary  condition  was  applied  to  FEM  model  as 
shown  below.  (Eq.  5) 

0{r,z  =  O,t)-2AD(d/dz)</>{r,z,t]^o  =0  (5) 


with 


A  =  (l  +  rrf)/(l-rrf) 


(6) 


where  rd  is  the  internal  reflectance,  due  to  refractive  index  mismatch  between  the  air  and  the  tissue.  This 
can  be  estimated  using  the  following  empirically  determined  equation. 


rd  =  -1.440/j"2  +0.710/T1  +  0.668  +  0.0636/j 


(7) 


with 


n  =  tissue  /  nair 


(8) 


We  used  1.4  and  1.0  as  values  for  ntissue  and  nair,  respectively,  and  then  we  could  get  3.25  as  a  value  for  A. 
We  also  used  0.033  cm  as  a  diffusion  coefficient  of  phantom  since  it  gives  reduced  scattering  coefficient, 
10  cm'1,  of  the  regular  tissue,  and  0.03  and  1.5  cm'1  were  used  as  absorption  coefficients  for  phantom  and 
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dye  in  VMDs,  respectively.  The  simulation  model  was  then  meshed  in  the  mesh  mode  resulted  in  1 147 
elements  and  609  nodes.  Finally,  it  is  solved  using  the  stationary  nonlinear  solver  type. 

Two  different  flow  rates  in  tumor  phantom  study  were  simulated  by  lengthening  the  heights  of 
Rl,  R2,  and  R3,  R4  with  different  frame  rates.  In  other  words,  a  fast  flow  rate  was  simulated  by 
increasing  the  heights  of  Rl  and  R2  with  a  large  increment  per  frame  (e.g.,  0  to  2.4  cm  with  0.6  cm 
increment  per  frame),  and  a  slow  flow  rate  was  simulated  by  increasing  heights  of  R3  and  R4  with  small 
increment  per  frame  (0  to  2.4  cm  with  0.02  cm  increment  per  frame).  Therefore,  there  existed  a  20  times 
difference  in  flow  rate  between  two  VMDs. 

We  also  changed  a  position  of  source  light  to  see  the  effect  of  vessel/VMD  geometry  on  bi-phasic 
feature  of  the  tumor  hemodynamic  measurements.  We  changed  the  position  of  light  source  to  the  right 
side  of  the  phantom  (Fig.  10a)  and  to  the  center  of  the  phantom  (Fig.  10b).  In  the  geometry  of  Fig.  10a, 
the  light  penetrates  the  slow-flow  VMD  first  and  then  the  fast-flow  VMD.  This  simulation  was  done  to 
find  whether  in  this  geometry,  the  light  signal  changes  still  have  the  double  exponential  behavior.  In  other 
words,  we  wish  to  investigate  whether  the  bi-phasic  hemodynamic  feature  depends  on  the  orientations  of 
the  vessels/VMD  networks. 

Multi-Channel  NIR  Spectroscopy 

Our  multi-channel  continuous  wave  NIR  spectroscopy  system  has  one  light  source  and  eight 
detectors,  but  due  to  the  finite  tumor  dynamic  phantom  size  and  for  the  simplicity  of  this  initial 
experiment,  we  used  only  three  detectors  to  monitor  light  absorption  changes  in  tumor  phantoms  during 
tumor  phantom  experiments.  Based  on  the  modified  Beer-Lambert’s  law,  the  data  presented  in  this  paper 
were  analyzed  using  amplitude  values  to  find  the  changes  in  absorption  (Eq.  9). 

Apa  =  Pa  B  -  Ha  T  =  log  (A0/AT)  /  L,  (9) 
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where  L  is  the  optical  path  length,  and  AB  and  ATare  baseline  and  transient  amplitudes  of  the  measured 
optical  signals,  respectively. 


RESULTS 

Tumor  Dynamic  Phantom  Study 

Figure  6a  shows  the  absorption  changes  (Apa,  wavelength=730nm)  of  tumor  phantom  at  three 
different  positions  when  dye  was  injected  into  VMDs.  The  first  step  for  this  experiment  was  injecting  dye 
into  only  VMD-2,  which  has  large  diameter  tubing  with  a  20  ml/hr  flow  rate.  As  we  can  see  from  the 
result,  the  reading  from  D 1  showed  largest  increase  of  Apa,  and  D2  showed  smallest  increase  of  Apa  since 
at  the  position  of  D 1 ,  it  would  detect  optical  signal  changes  mostly  from  the  VMD-2.  The  Apa  increase  of 
D3  was  between  those  from  D1  and  D2,  as  we  expected.  A  similar  pattern  is  also  shown  at  step  2  where 
there  was  a  dye  flow  (20  ml/hr)  only  into  VMD-1,  which  has  small  diameter  tubing.  At  this  step,  D2  had 
a  largest  increase  of  A|i.a  and  D1  showed  smallest  Apa  increase  with  the  same  reason  that  was  mentioned 
above.  Step  1  and  2  clearly  showed  us  that  the  D1  and  D2  are  placed  on  the  tumor  phantom  properly  for 
this  experiment  as  we  designed,  and  clearly  demonstrates  that  signal  from  heterogeneity  can  be  detected 
stronger  if  it  is  near  the  detector.  The  third  step  in  the  measurement  was  to  inject  dye  into  both  VMD-1 
and  VMD-2  simultaneously  with  a  same  flow  rate  as  before  (20  ml/hr).  At  this  step,  we  can  observe  two 
features:  1)  the  time  profile  taken  at  D2  has  a  faster  transition  time  than  that  at  Dl,  and  2)  the  time  profile 
taken  at  D3  has  a  bi-exponential  characteristic,  which  is  very  similar  to  those  we  often  observed  in  our 
breast  tumor  dynamic  measurements.25 

To  understand  these  phantom  measurements,  let  us  consider  the  velocities  at  each  tube.  Even 
though  the  flow  rate  was  the  same  for  both  of  the  tubes,  the  velocities  of  the  ink  solutions  in  the  two 
VMDs  networks  were  different  due  to  the  different  diameters.  The  solution  in  the  small  diameter  tube 
had  a  velocity  2.74  times  faster  than  that  of  the  large  diameter  tube,  based  on  the  following  relationship: 
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VVMP- 1  _  rVMD-2 
VVMD- 2  rVMD- 1 


0.42 

0.2552 


2.46 


(10) 


where  vvMD/.and  vVmd-2  are  velocities  of  the  ink  solution  in  the  VMD-1  (small  diameter  tubing)  and 
VMD-2  (large  diameter  tubing),  respectively,  and  rVMD  i,  and  rVMD-2  present  the  radii  of  the  VMD-1  and 
VMD-2,  respectively.  This  velocity  difference  may  be  the  reason  why  the  transition  time  is  faster  in  the 
VMD-1  than  that  in  the  VMD-2.  When  we  increased  the  flow  rate  of  large  diameter  tube,  we  could  not 
see  the  two  exponential  behaviors  from  D3.  Furthermore,  since  D3  was  in  the  transmission  geometry  and 
was  located  near  both  of  the  VMDs,  the  signal  obtained  at  D3  should  reflect  the  dynamic  changes 
occurred  in  both  of  the  VMDs.  We  expect  that  the  bi-exponential  feature  recorded  by  D3  results  from  a 
superposition  between  the  two  different  dynamic  transitions  at  the  VMD-1  and  VMD-2. 

For  the  further  analysis,  Ajxa  increases  detected  from  three  detectors  at  step  3  in  Figure  6a  was 
fitted  by  using  our  mathematical  models  to  obtain  amplitudes  and  time  constants  values.  (Figure  6b)  The 
fitting  values  for  each  curve  are  shown  in  the  table  1.  This  result  shows  that  when  dye  flows  into  both 
VMD-1  and  VMD-2,  fast  time  constant  (t1=1  1.71  ±  3.49)  and  slow  time  constant  (t2=35.17  ±  7.25)  from 
D3  are  very  close  to  xl  (=7.3  ±  0.25)  at  D2  (VMD-1)  and  tl  (=42.8  ±  1.53)  at  D1  (VMD-2),  respectively. 
This  indeed  shows  that  two  exponential  behaviors  of  tumor  blood  oxygenation  during  carbogen  inhalation 
can  be  from  two  different  blood  vessels  diameters  in  tumors. 

The  phantom  1  data  presented  here  basically  demonstrate  that  we  can  experimentally  mimic  the 
bi-exponential  behavior  with  two  different  velocities,  or  perfusion  rates,  from  two  different  diameter 
VMDs.  This  knowledge  may  help  us  understand  the  bi-phasic  features  observed  in  the  breast  tumor 
vascular  oxygenation  measurements  under  carbogen  inhalation.  Furthermore,  the  next  phantom  2 
experimental  results  demonstrate  that  the  bi-exponential  dynamics  may  also  result  from  different  flow 
rates,  given  the  same  diameters  for  the  two  VMDs. 

Figure  7a  shows  the  experimental  results  from  tumor  phantom  2,  and  this  phantom  is  different 
from  tumor  phantom  1  by  having  two  VMD-ls  instead  of  one  VMD-1  and  one  VMD-2  in  tumor  phantom 
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1.  Since  we  now  have  two  VMD-ls,  the  velocities  of  dye  in  VMDs  will  be  totally  depending  on  the  dye 
flow  rates  controlled  by  two  syringe  pumps.  The  source  and  detector  positions  were  same  as  Figure  6a, 
which  is  shown  in  the  Figure  4b,  and  we  observed  the  results  similar  to  those  shown  in  Figure  6a  by 
alternating  flow  rates. 

The  first  step  for  this  experiment  was  injecting  dye  only  into  bottom  VMD-1  with  a  flow  rate  of 
20  ml/hr.  As  a  result,  D2  showed  largest  increase  of  Apa  while  D1  showed  smallest  increase  of  A|a.a 
ecause  D1  was  located  quite  away  from  the  bottom  VMD-1.  For  the  second  step,  top  VMD-1  was 
injected  with  dye  at  the  same  flow  rate  of  20  ml/hr.  Here,  D 1  showed  largest  increase  of  Apa  and  D2 
showed  smallest  increase  of  Apa.  In  both  of  the  above  cases,  the  signals  from  D3  showed  similar  profiles 
as  the  others,  without  clear  bi-exponential  patterns.  Then,  as  the  third  step,  in  dye  was  injected  into  both 
top  and  bottom  VMD-ls  with  the  same  flow  rate  (20  ml/hr).  Similarly,  we  did  not  observe  any  clear 
features  of  the  two-exponential  increase  of  Apa.  For  the  fourth  step,  wee  injected  dye  into  both  VMD-ls, 
but  with  two  different  flow  rates:  5  ml/hr  for  top  VMD-1  and  20  ml/hr  for  the  bottom  VMD-1.  Now  we 
can  see  bi-exponential  behavior  appears  on  D3  since  D3  is  detecting  signals  from  both  top  and  bottom 
VMD-ls,  which  have  two  different  dye  flow  rates. 

Once  again,  time  constant  analysis  was  performed  for  the  increase  of  Apa  from  three  detectors  at 
step  4  in  Fig.  7a  and  is  shown  in  Fig.  7b.  This  time  constant  analysis  shows  that  Apa  increases  from  D1 
and  D2  are  well  fitted  by  mono  exponential  model  while  Apa  increases  from  D3  is  fitted  well  with  double 
exponential  models.  As  we  can  see  from  the  table  2,  which  the  values  of  fitted  parameters  are 
summarized,  fast  time  constant  (xl=  18.13  ±  0.87)  and  slow  time  constant  (x2=  133.66  ±  1 1.36)  from  D3 
are  well  matched  with  xl  in  D2  (20.92  ±  0.53)  and  xl  in  D1  (131.19  ±  3.46).  This  suggests  that  the 
dynamic  signals  obtained  from  D3  with  a  fast  and  slow  component  indeed  result  from  the  dynamic 
changes  in  both  of  the  VMD-ls. 

FEM  Simulation 
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Figure  8  shows  an  example  of  the  FEM  tumor  phantom  simulation  to  further  support  our  tumor 
phantom  2  experiment.  Figure  8a  is  the  light  distribution  of  dynamic  phantom  where  there  was  dye 
injection  into  only  R1  and  R2  (bottom  VMD-1  in  Fig.  7a)  with  a  fast  flow  rate  and  Fig.  8b  is  the 
simulation  when  dye  flowed  only  into  R3  and  R4  (top  VMD-1  in  Fig.  7a)  with  a  slow  flow  rate.  Finally, 
Figure  8c  shows  how  light  distributes  in  a  tumor  dynamic  phantom  when  dye  was  injected  into  both 
VMD-ls  with  two  different  flow  rates.  As  we  mentioned  earlier  at  materials  and  methods  section, 
different  flow  rate  was  mimicked  by  increasing  the  length  of  light  absorbing  object  (Rl,  R2,  R3  and  R4) 
with  the  different  frame  rates.  Figure  8a  is  the  result  of  simulation  at  6th  frame  while  Fig.  8b  is  the 
simulated  result  at  120th  frame.  This  implies  that  Rl  and  R2  have  20  times  faster  flow  rate  than  that  in  R3 
and  R4.  Figure  8c  is  also  the  result  of  simulation  at  120th  frame  when  dye  flows  in  both  VMD-ls. 

For  comparison  with  the  results  from  phantom  experiments,  we  extracted  light  intensity 
values  at  three  positions  (2,0),  (-2,0),  and  (0,2)  from  each  frame  of  simulations  to  calculate  Apa,  which  are 
shown  in  the  figure  9.  These  three  positions,  (2,0),  (-2,0),  and  (0,2),  are  corresponding  to  Dl,  D2  and  D3 
in  Fig.  7a,  and  time  in  these  plots  were  obtained  by  multiplying  2  seconds  to  each  frame.  Therefore,  Fig. 
8a  will  have  12  seconds  to  reach  the  maximum  Apa  since  it  has  only  6  frames  to  simulate  fast  flow,  and 
Fig.  8b  and  8c  will  have  240  seconds  to  have  a  maximum  absorption  because  they  have  120  frames  of 
slow  flow  simulation. 

Figure  9a  and  9b  show  Apa  from  three  positions  of  phantom  in  fast  flow  only  simulation  and  slow 
flow  only  simulation,  respectively.  Figure  9a  shows  the  Ajia  is  largest  in  (-2,0)  position  and  smallest  in 
(2,0)  position  when  there  is  a  dye  flow  only  into  fast  flow  VMD-1,  and  this  result  is  matching  with  the 
results  from  Fig.  6a  and  7a.  The  similar  results  could  be  obtained  when  dye  flowed  only  into  slow  VMD- 
1  as  we  can  see  from  the  Fig.  9b.  Figure  9c  shows  Apa  from  three  positions  in  Figure  8c,  and  the  Ap„  at 
(0,2)  position  shows  bi-phasic  behavior  as  we  could  see  from  animal  tumors  and  our  tumor  dynamic 
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phantoms.  This  indeed  supports  that  it  is  necessary  to  have  a  different  flow  rate  for  tumor  to  have  tri¬ 
phasic  blood  oxygenation  dynamics  during  carbogen  inhalations. 

We  also  changed  the  position  of  source  light  to  the  right  side  of  the  phantom  (Fig.  10a)  and  center 
of  the  phantom  (Fig.  10b)  to  see  how  blood  vessel/VMD  geometry  in  tumor/phantom  affects  the  bi-phasic 
feature  of  the  tumor  hemodynamics  during  carbogen  inhalation.  As  described  in  Figure  5,  R1  and  R2 
represent  the  fast  flow  VMD  while  R3  and  R4  represent  the  slow  flow  VMD.  Light  intensity  values  at 
several  positions  were  also  extracted  to  calculate  A(i.,  and  were  plotted  in  Figure  1  la  and  lib.  Figure  1  la 
shows  that  (0,-2)  and  (0,2)  positions  detect  little  signals  from  fast  flow  VMD  since  both  positions  will 
mostly  detect  light  from  the  slow  flow  VMD  with  this  geometry.  However,  this  geometry  allows  (-2,0) 
position  can  detect  signals  from  both  slow  and  fast  flow  VMDs  so  that  we  can  see  bi-phasic  behavior 
from  (-2,0)  position.  Figure  1  lb  also  plotted  the  A|ia  obtained  from  four  positions  on  tumor  phantom 
simulation  (Fig.  10b)  and  showed  that  (-2,0)  and  (2,0)  positions  were  getting  signals  mostly  from  only 
fast  flow  and  slow  flow  VMD,  respectively.  However,  (0,2)  and  (0,-2)  positions  were  showing  bi-phasic 
behavior  because  they  were  detecting  light  signals  from  both  fast  flow  and  slow  flow  VMDs.  All  of  these 
results  prove  that  there  will  be  bi-phasic  features  of  tumor  hemodynamics  as  long  as  two  different 
perfusion  rates  exist  in  tumors  regardless  of  vessel  geometry. 

DISCUSSION 

Both  tumor  dynamic  phantom  experiments  and  FEM  simulation  results  supported  our  hypothesis 
that  the  bi-phasic  tumor  hemodynamic  feature  results  from  a  well-perfused  and  poorly  perfused  region  in 
the  tumor  vasculature  under  carbogen  inhalation.  Through  our  phantom  study,  we  were  able  to  find  that 
main  cause  for  having  a  bi-phasic  behavior  is  coming  from  different  blood  flow  rate,  more  precisely  from 
different  blood  velocity  in  tumor.  These  differences  in  velocity  can  be  from  different  blood  vessel 
diameters  with  same  blood  flow  rate  or  from  different  blood  flow  rate  with  same  vessel  diameter,  which 
tumor  vasculatures  have  both.  Tumor  blood  vessels  are  known  to  be  very  leaky,  longer  in  vessel  lengths, 
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larger  in  vessel  diameter,  and  also  their  local  microvessel  density  is  much  more  heterogeneous  compare  to 
normal  tissues.8  In  addition,  solid  tumors  usually  develop  hypoxia,  which  can  be  from  poor  perfusion  in 
the  central  region  when  they  grow  bigger.  Meanwhile,  the  peripheral  region  of  tumor  is  normally  well 
perfused  so  that  it  can  be  well  provided  with  nutrition  and  oxygen.  Therefore,  tumor  vasculature  can  be  a 
mixed  structure  between  well-perfused  region  and  poorly  perfused  region,  which  can  be  observed  as  a  bi- 
phasic  feature  of  hemodynamics  with  respiratory  challenges.  Therefore,  the  multi-channel  NIRS  system 
will  allow  us  to  study  the  dynamic  heterogeneity  at  different  locations,  showing  that  the  tumor 
vasculature  in  the  poorly  perfused  region  (with  a  lower  perfusion  rate)  can  be  eventually  oxygenated  if  the 
two  dynamic  components  are  observed.  On  the  other  hand,  if  only  the  fast  component  exists,  it  suggests 
that  the  poorly  perfused  region  in  the  tumor  may  not  have  much  chance  to  be  improved  for  its  vascular 
oxygenation  under  carbogen  intervention. 

FEM  simulation  results  strongly  supported  our  tumor  dynamic  phantom  experimental  data  and 
showed  that  the  bi-phasic  behavior  always  existed  regardless  of  geometry  of  light  source  and  VMDs  in 
phantom  as  long  as  a  detector  is  located  somewhere  it  can  collect  signals  from  both  fast  flow  and  slow 
flow  VMDs.  To  compare  the  results  from  tumor  phantom  experiments  with  FEM  simulation  results,  we 
divided  Apa  by  maximum  value  of  Apa  from  each  detector  for  normalization.  (Figure  12a  and  12b) 
Figure  12a  is  a  normalized  plot  of  step  4  in  Fig.  7a  (tumor  phantom  experiment)  and  Fig.  12b  is  another 
normalized  plot  of  Fig.  9c.  (FEM  simulation)  These  figures  show  how  much  FEM  simulation  results  are 
close  to  tumor  phantom  experiments.  Both  of  them  show  that  Apa  increases  from  D3  and  (0,2)  position 
will  be  well  fitted  by  using  a  double  exponential  model.  The  early  time  of  Apa  increases  from  D3  and 
(0,2)  position  (inside  the  circle)  have  larger  Apas  than  D2,  Dl,  (-2,0)  and  (2,0),  which  can  be  explained  by 
locations  of  detectors.  Since  the  dye  flow  started  from  D3  or  (0,2)  side,  both  D3  and  (0,2)  position  could 
detect  Apa  changes  earlier  and  also  larger  than  those  from  the  other  detectors  (Dl,  D2,  (2,0),  and  (-2,0)), 
which  were  located  a  little  far  from  the  dye  flow  starting  position.  These  results  indeed  show  that  FEM 


32 


simulation  can  well  explain  the  results  from  tumor  phantom  experiments  by  using  light  diffusion 
equation. 

As  we  described  in  materials  and  methods  section  earlier,  Y1/Y2  may  be  associated  with  vascular 
volume  or  density  of  two  regions,  and  f\lf2  is  related  to  the  ratio  of  blood  perfusion  rates  between  region  1 
and  2.  Our  tumor  phantom  experiments  and  FEM  simulation  of  tumor  dynamic  phantom  showed  that  A 
or  y  values  are  related  with  absorption  changes  of  VMDs  in  tumor  phantoms,  and  T  or/are  related  to  the 
velocity  of  dye  flow.  The  intensity  of  absorption  changes  that  measured  from  detectors  depends  on  the 
absorption  coefficient  of  dye,  the  tube  length  or  diameter  in  VMD,  wrapping  number  of  small  tubing 
around  big  tubing,  and  location  of  detectors.  Time  constant  and  perfusion  rate  are  mainly  depending  on 
dye  flow  rate,  but  it  can  be  also  affected  by  diameter  and  length  of  tube  in  VMD.  In  tumor’s  case,  pa  of 
dye  corresponds  to  the  level  of  blood  oxygenation,  and  tube  length  or  diameter  can  be  thought  as  those  of 
blood  vessel.  The  number  of  winding  small  tubing  around  big  tubing  can  be  thought  as  a  vascular  density 
in  tumor,  which  will  affect  the  values  of  both  y  and  /. 

Our  previous  measurements  of  animal  tumors  were  performed  in  a  transmittance  mode  by  one 
channel  NIRS,  which  could  obtain  only  global  measurements  of  tumor  hemodynamics.  Therefore,  those 
results  could  not  show  intratumoral  or  intertumoral  vascular  heterogeneities.  However,  multi-channel 
NIRS,  by  comparing  Y1/Y2  and  fiff2  among  the  signals  taken  at  different  locations  from  the  same  tumor  or 
from  different  tumors,  will  enable  us  to  study  intratumoral  or  intertumoral  heterogeneity.  Since  tumor 
vascular  structure  will  be  modified  during  therapies,  we  will  be  able  to  observe  the  changes  of  Y1/Y2  and 
f\/f2  through  time  by  giving  respiratory  challenges  such  as  carbogen  or  pure  oxygen  inhalation.  Therefore, 
multi-channel  NIRS  or  imaging  not  only  can  detect  vascular  heterogeneity  of  tumor  but  also  can  be  a 
prognostic  tool  for  monitoring  the  early  effects  of  tumor  treatments  such  as  radiotherapy,  photodynamic 
therapy,  or  chemotherapy.  The  detection  of  early  effects  of  cancer  therapies  will  benefit  patients  to  have 
a  higher  survival  rate  by  treating  the  patients  with  proper  therapies. 


33 


CONCLUSION 


We  developed  tumor  phantom  models  to  demonstrate  our  hypothesis  on  tumor  hemodynamics  during 
carbogen  inhalation.  We  believe  that  bi-phasic  feature  of  tumor  blood  oxygenation  increases  during 
carbogen  inhalation  is  from  tumors’  distinct  vascular  structure,  which  is  composed  of  well-perfused 
region  and  poorly  perfused  region.  Our  tumor  phantom  experiments  were  performed  to  find  what  can 
cause  tumor  hemodynamics  to  have  two  time  constants  when  it  is  fitted  with  our  mathematical  model. 
We  found  that  two  time  constants  in  tumor  hemodynamic  models  can  be  from  different  blood  flow 
velocities  or  anything  that  can  cause  changes  in  blood  flow  velocities  such  as  blood  vessel  diameter,  and 
length.  In  addition,  we  also  found  that  two  different  amplitudes  in  tumor  hemodynamic  models  are  from 
differences  in  absorption  in  two  region,  which  are  possibly  due  to  different  blood  oxygenation  level 
(oxyhemoglobin  concentration),  or  vascular  density.  FEM  simulation  using  a  light  diffusion  equation  was 
performed  to  simulate  our  tumor  phantom  experiments.  These  simulated  results  were  well  matched  with 
our  tumor  phantom  experiments,  which  indeed  supported  our  hypothesis. 

NIRS  is  a  portable,  low  cost,  and  real  time  measurement  system  that  can  monitor  changes  of 
vascular  oxygen  levels  by  using  two  wavelengths.  We  have  previously  used  a  single-channel  NIRS 
system  with  one  light  source  and  one  detector  for  global  measurements  of  A[Hb02]  and  A[Hbtotai]  in 
tumors  during  respiratory  challenges.25  Now,  this  study  demonstrates  that  the  NIRS  multi-channel 
approach  has  the  great  potential  to  detect  and  monitor  tumor  heterogeneity  under  therapeutic  or 
adjuvant  interventions.  In  the  future  we  propose  1)  to  further  investigate  and  understand  the  meaning 
of  vasculature  coefficient,  y,  and  2)  to  develop  an  NIR  imaging  system  to  image  tumor  oxygenation 
heterogeneities  under  various  interventions  used  to  improve  the  efficacy  of  cancer  treatments. 
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Table  1.  Summary  of  fitted  parameters  obtained  at  the  three  detectors  from  step  3  in  Figure  6a 


Parameters 

Mono-Exponential  fitting 

AHb02  =  A,[l-exp(-t/Ti)] 

Double-Exponential  fitting 

AHb02  =  A1[l-exp(-t/T,)]  + 

A2[l-exp(-t/x2)] 

Detectors 

Detector  D1 

Detector  D2 

Detector  D3 

A,  (mM/DPF) 

0.49  ±0.007 

0.34  ±0.001 

0.079  ±0.041 

Ti  (min) 

42.8+1.53 

7.3  ±  0.25 

11.71  ±3.49 

A2  (mM/DPF) 

0.14  ±0.038 

x2  (min) 

35.17  ±7.25 

5 

X 

0.039 

0.016 

0.0018 

R 

0.99 

0.98 

0.99 
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Table  2.  Summary  of  fitted  parameters  obtained  at  the  three  detectors  from  step  4  in  Figure  7a 


Parameters 

Mono-Exponential  fitting 

AHb02  =  A[[l-exp(-t/Xi)] 

Double-Exponential  fitting 

AHb02  =  A|[l-exp(-t/T,)] + 

A2[l-exp(-t/x2)] 

Detectors 

Detector  D1 

Detector  D2 

Detector  D3 

A,  (mM/DPF) 

0.26  ±  0.003 

0.1710.001 

0.068  ±  0.0024 

ti  (min) 

131.19  ±3.46 

20.92  ±  0.53 

18.1310.87 

A2  (mM/DPF) 

0.047  1  0.0018 

x2  (min) 

133.661  11.36 

r 

0.057 

0.019 

0.0018 

R 

0.98 

0.97 

0.99 
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FIGURE  CAPTIONS 


Figure  1. Results  obtained  with  the  NIR  instrument  from  a  4.5-cm3  rat  breast  tumor  while  the  breathing 
gas  was  switched  from  33%  O2  to  carbogen.  The  thicker  solid  curve  represents  AHbCF,  the  thinner  solid 
curve  represents  AHb„  and  the  dashed  curve  with  the  filled  circles  represents  arterial  saturation.  (Liu  et  al. 
2000)25 

Figure  2.  Schematic  diagram  of  light  transmitting  patterns  in  tumor  when  tumor  has  two  distinct  perfusion 
regions.  Center  of  tumor  represents  poorly  perfused  region  and  peripheral  region  of  tumor  with  gray  color 
is  representing  a  well-perfused  region.  Since  the  light  transmitting  volume  is  different  in  each  detector, 
each  detector  will  show  the  different  fitted  parameters. 

Figure  3.  Schematic  diagram  for  one  embedded  vascular  modeling  device  (VMD):  the  small  tubing 
wound  outside  of  core  tubing  simulates  blood  vessels  within  a  breast  tumor. 

Figure  4.  Experimental  setup  for  the  tumor  dynamic  phantom  study,  (a)  Two  syringe  pumps  were 
connected  to  two  VMDs  in  tumor  dynamic  phantom  for  controlling  dye  flow  rates.  Light  was  transmitted 
from  light  source  through  tumor  phantom  and  collected  at  three  different  detectors  for  data  processing  in 
computer  (b)  Enlarged  tumor  dynamic  phantom  embedded  with  two  VMDs.  Tumor  phantom  1  has  VMD- 
1  and  VMD-2  as  shown  here  and  tumor  phantom  2  has  two  VMD- Is. 

Figure  5.  The  geometry  for  FEM  simulations  of  a  tumor  dynamic  phantom  2.  R1  and  R2  rectangles 
represent  fast  flow  VMD  while  R3  and  R4  represent  slow  flow  VMD. 

Figure  6.  (a)  3-channel  NIRS  results  measured  on  a  tumor  dynamic  phantom  1  with  two  size  VMDs 
(VMD-1  and  VMD-2).  Three  traces  represent  the  readings  at  D3  (in  transmission  mode  and  located 
between  the  two  VMDs),  D1  (near  the  VMD-2),  and  D2  (near  the  VMD-1).  (b)  Temporal  profiles  of  the 
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NIRS  measurements  from  Dl,  D2,  and  D3  with  one-exponential  fitting  (Dl,  and  D2)  and  bi-exponential 
fitting  (D3). 

Figure  7.  (a)  3-channel  NIRS  results  measured  on  a  tumor  dynamic  phantom  2  with  two  VMD-ls.  Three 
traces  represent  the  readings  at  D3  (in  transmission  mode  and  located  between  the  two  VMDs),  Dl  (near 
the  top  VMD-1),  and  D2  (near  the  bottom  VMD-1).  (b)  Temporal  profiles  of  the  NIRS  measurements 
from  Dl,  D2,  and  D3  with  one-exponential  fitting  curves  (Dl,  and  D2)  and  bi-exponential  fitting  curve 
(D3). 

Figure  8.  FEM  simulations  of  light  distribution  inside  of  a  dynamic  phantom,  (a)  The  results  with  a  fast 
simulated  flow  rate  only,  (b)  the  results  with  a  only  slow  flow  rate,  and  (c)  the  results  with  both  fast  and 
slow  simulated  flow  rates. 

Figure  9.  Absorption  changes  measured  at  three  locations  in  FEM  simulation  (Fig.  8)  during  fast  flow 
only  (a),  slow  flow  only  (b),  and  both  fast  and  slow  flow  (c). 

d 

Figure  10.  FEM  simulations  of  light  distribution  inside  of  a  dynamic  phantom,  (a)  The  results  when  the 
light  source  is  located  on  (2,0)  position,  (b)  the  results  when  the  light  source  is  located  in  the  center  of  the 
phantom. 

Figure  11.  Absorption  changes  measured  at  three  locations  in  FEM  simulation  (Fig.  10a)  during  both  fast 
and  slow  flow  (a),  and  Apa  measured  at  four  locations  in  FEM  simulation  (Fig.  10b)  during  both  fast  and 
slow  flow  (b). 
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Figure  12.  Normalized  Ap.,  plots  are  shown  to  compare  the  results  from  between  tumor  dynamic  phantom 
experiment  and  FEM  simulation:  (a)  Normalized  A|ia  plots  of  step  4  in  Fig.  7a,  (b)  Normalized  Apa  plots 
of  Fig.  9c. 
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ABSTRACT 

In  our  previous  report,  we  showed  that  there  was  a  severe  heterogeneity  of  tumor  tissue  oxygen  dynamics 
even  in  a  same  tumor  by  using  l9F  MR  p02  mapping  and  near-infrared  spectroscopy  (NIRS).  In  this  study,  we 
applied  a  three  channel  NIRS  system  to  reveal  heterogeneity  of  tumor  vascular  oxygen  dynamics  during 
respiratory  challenges.  When  the  implanted  tumor's  size  reached  3  cm,  three  photo  detectors  were  attached  on 
the  surface  of  tumors  in  three  different  positions  to  measure  changes  in  oxygenated  hemoglobin  concentration 
([HbOi])  while  the  inhaled  gas  was  alternated  between  air  and  carbogen  (95%  oxygen,  5%  C02).  Significant 
changes  in  tumor  oxygenation  were  observed  accompanying  respiratory  challenge,  and  these  changes  could 
be  modeled  with  either  one  or  two  exponential  components  with  fast  and  slow  time  constants.  Time  constant, 
amplitude,  vascular  coefficient,  and  the  ratio  of  perfusion  rate  obtained  from  fitting  curves  of  A[Hb02] 
measured  from  three  photo  detectors  revealed  that  there  were  indeed  different  responses  of  tumor  vascular 
oxygenation  during  carbogen  inhalation  from  three  different  locations.  These  results  clearly  show  that  tumors 
are  highly  heterogeneous  in  terms  of  vascular  oxygenation  response  to  carbogen  inhalation.  This  study 
demonstrates  that  the  NIR  technology  can  provide  an  efficient,  real-time,  non-invasive  way  for  monitoring 
tumor  physiology,  and  it  may  have  prognostic  value  and  promised  insight  into  tumor  vascular  development 
and  angiogenesis. 

Keywords:  Tumor  vascular  oxygenation,  tumor  heterogeneity,  NIR  Spectroscopy,  carbogen 
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1.  INTRODUCTION 

Solid  tumors  develop  regions  of  hypoxia  during  their  growth  due  to  an  imbalance  between  the  rate  of  tumor 
cell  proliferation  and  the  proliferation  and  branching  of  the  blood  vessels.1'2  3  Tumor  hypoxia  is  responsible 
for  the  failure  of  radiotherapy,4  some  forms  of  chemotherapy,5  and  photodynamic  therapy.6  In  addition,  a 
number  of  clinical  trials  have  found  that  patient  survival  measured  either  as  tumor  regression  or  as  local 
control  depends  largely  on  tumor  oxygenation.7  Tumor  hypoxia  can  occurr  through  diffusion-limited  or 
chronic  hypoxia  and  perfusion-limited  or  acute  hypoxia.8  Therefore,  measurement  of  distribution  of  tumor 
perfusion  rate  could  be  important  for  tumor  treatment  planning  and  assessing  methods  designed  to  modulate 
tumor  oxygenation. 

Solid  tumors  are  known  to  exhibit  heterogeneous  blood  flow  distribution.9'10  There  are  many  methods  to 
measure  tumor  perfusion  heterogeneity,  such  as  Doppler  ultrasound,11  dynamic  contrast  MRI,12  and  the  use 
of  tumors  grown  in  windowed  chambers.13  Intensive  studies  from  Mason’s  group  using  19F  MR  p02  mapping 
have  revealed  intratumoral  heterogeneity  of  p02  distribution  and  also  heterogeneous  response  to  hyperoxic 
gas  breathing.14'15'16  The  group’s  findings  of  severe  p02  heterogeneity  in  tumors  can  indirectly  indicate  the 
heterogeneous  distribution  of  blood  flow,  since  tissue  p02  level  is  decided  by  a  balance  between  the  supply  of 
oxygen  from  blood  vessels  and  oxygen  consumption  by  tissue  cells.17  Unlike  the  19F  MR  p02  mapping 
technique,  near-infrared  spectroscopy  (NIRS)  techniques  cannot  measure  p02  in  tissue  but  can  measure  in 
vivo  hemoglobin  saturation  and  concentration.  NIRS  has  been  applied  to  muscle,18 1920  brain,21 22  23  and 

24,25,26 

tumor.  ’ 

In  our  previous  study,  we  established  a  mathematical  model  of  tumor  oxygen  dynamics  during  hyperoxic  gas, 
such  as  carbogen  (95%  C02  and  5%  02)  or  oxygen  inhalation.26  For  our  model,  we  formed  a  hypothesis  that 
changes  in  oxygenated  hemoglobin  concentration  (A[Hb02])  signal  is  composed  of  signals  from  well- 
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perfused  and  poorly  perfused  region  to  explain  why  there  are  two  different  time  constants  to  fit  the  A[Hb02] 
data  using  our  model.  In  addition  to  time  constant,  we  could  get  amplitude  of  fitted  curve,  ratio  of  vascular 
coefficients,  and  ratio  of  perfusion  rates  by  applying  our  dynamic  tumor  oxygenation  model.  Therefore,  we 
will  be  able  to  detect  intratumoral  perfusion  heterogeneity  by  applying  multi-channel  NIRS,  since  each 
detector’s  signal  will  come  from  a  different  region  of  tumor  and,  thus,  it  will  show  differences  of  curve  fitting 
parameters  among  multi-detectors. 

The  goals  of  this  study  were  to  measure  A[Hb02]  in  tumor  vasculature  during  carbogen  inhalation  using  a 
three  channel  NIRS  system,  and  to  reveal  intratumoral  vascular  heterogeneity  by  fitting  the  A[Hb02]  data 
with  our  mathematical  model. 

2.  MATERIALS  AND  METHODS 

2.1  Tumor  Model  and  Measurement 

For  our  study,  murine  mammary  adenocarcinomas  13762NF  grown  in  the  hind  limb  of  adult  female  Fisher 
344  rats  (-200  g)  were  used.  Once  these  tumors  reached  approximately  2-3  cm  in  diameter,  the  rats  were 
anesthetized  with  0.2  ml  ketamine  hydrochloride  (100  mg/ml;  Aveco,  Fort  Dodge,  IA)  and  maintained  under 
general  gaseous  anesthesia  using  a  small  animal  anesthesia  unit  with  air  (1  dm3/ min)  and  1.3%  isoflurane 
(Ohmeda  PPD  Inc.,  Fort  Dodge,  IA)  through  a  mask  placed  over  the  mouth  and  nose.  After  anesthesia,  the  rat 
was  placed  on  a  warm  blanket  to  maintain  body  temperature,  which  was  monitored  with  a  rectally-inserted 
thermal  probe  connected  to  a  digital  thermometer  (Digi -Sense,  model  91100-50,  Cole-Parmer  Instrument 
Company,  Vernon  Hills,  IL).  Tumors  were  shaved  to  improve  optical  contact  for  transmitting  light,  and  a 
light  source  and  three  detectors  were  attached  to  the  tumor  using  posts  and  swivel  post  clamps.  A  pulse 
oximeter  (model:  8600V,  Nonin,  Inc.)  was  placed  on  the  hind  foot  to  monitor  arterial  oxygenation  (Sa02)  and 
heart  rate.  Figure  1  shows  schematic  setup  for  animal  experiment  using  multi  channel  NIRS  system. 
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All  measurements  were  performed  in  the  dark  room,  and  measurements  were  initiated  while  the  rats  breathed 
air  for  10  minutes  to  get  a  stable  baseline.  After  10  min  of  baseline  measurement,  the  inhaled  gas  was 
switched  to  carbogen  for  15  minutes,  and  then  back  to  air  for  15  minutes.  Tumor  volume  V  (in  cm')  was 
estimated  as  V  =  (4it/3)  [(L+W+H)/6],27  where  L,  W,  and  H  are  the  three  respective  orthogonal  dimensions. 
Raw  amplitude  data  from  three  detectors  were  recorded  simultaneously  during  the  experiments,  and 
processed  after  the  experiments  to  obtain  the  changes  in  oxygenated  hemoglobin  concentration  (A[HbC>2])  and 
total  hemoglobin  concentration  (A[Hb,ota|]).  Time  constant,  amplitude,  Y1/Y2,  and /1//2  were  calculated  by  fitting 
the  data  using  Kaleidagraph  software  (Synergy  Software,  Reading,  PA). 

2.2  Multi-Channel  NIR  Spectroscopy 

Our  multi-channel  NIR  spectroscopy  system  has  one  light  source  and  eight  detectors  to  measure  light  signals 
from  eight  different  locations,  but  due  to  tumor  versus  detector  size  constraints  and  to  simplify  the 
experiment,  we  only  used  three  detectors  to  monitor  tumor  vascular  oxygen  dynamics  during  respiratory 
challenges..  In  common  with  our  previous  work26'28,  we  assume  that  oxyhemoglobin  and  deoxyhemoglobin 
are  the  only  significant  absorbance  materials  in  the  blood.  The  absorption  coefficients  comprise  the  extinction 
coefficients  for  deoxyhemoglobin  and  oxyhemoglobin  multiplied  by  their  respective  concentrations  (Eqs.  1 
and  2). 


=  eHb730[Hb]  +  eHbO2730[HbO2], 

(1) 

=  eHb850[Hb]  +  eHbO2850[HbO2]. 

(2) 

The  data  presented  in  this  paper  were  analyzed  using  modified  Beer-Lamberf  s  law  and  amplitude  values  to 
find  the  changes  in  absorption  (Eq.  3).  By  manipulating  Equations  1-3,  changes  in  oxygenated  hemoglobin. 
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deoxygenated  hemoglobin  and  total  hemoglobin  concentrations  were  calculated  from  the  transmitted 
amplitude  of  the  light  through  the  tumor  (Eqs.  4,  5  and  6). 

PaB  *  HaT  =  log  (AB/AT)  /  L, 

A|Hb02]  =  [-0.674*log  (AB/AT)73°  +  1.117*log  (AB/AT)85°J  /  L, 

A[Hb]  =  [0.994*log  (AB/AT)73°  -  0.376*log  (AB/AT)85°]  /  L, 

730  850 

A[Hb]l0Ul  =  A[Hb]  +  A[HbO2]=[0.32*log  (AB/\)  +  0.74 1  *log  (Ab/Aj.)  ]  /  L, 

(6) 

where  AB  =  baseline  amplitude;  AT  =  transition  amplitude;  L  =  optical  pathlength  between  source/detector. 
The  constants  were  computed  with  the  extinction  coefficients  for  oxy  and  deoxyhemoglobin  at  the  two 
wavelengths  used.29  In  principle,  L  should  be  equal  to  the  source-detector  separation,  d,  multiplied  by  a 
differential  pathlength  factor  (DPF),  i.e.,  L=d*DPF.  Little  is  known  about  DPF  for  tumors,  although  a  DPF 
value  of  2.5  has  been  used  by  others.30  Since  our  focus  is  on  dynamic  changes  and  relative  values  of  tumor 
[Hb02]  with  respect  to  carbogen  intervention,  we  have  taken  the  approach  of  including  the  DPF  in  the  unit, 
i.e.,  modifying  eq.  (4)  as  follows: 

A[Hb02]  =  [-0.674*log  (AB/AT)73°  +1.1 17*Iog  (Ag/A/50]  /  d,  (7) 

where  d  is  the  direct  source-detector  separation  in  cm,  and  the  unit  of  A[Hb02]  in  Eq.  (7)  is  mM/DPF. 

2.3  Mathematical  Model  of  Tumor  Vascular  Oxygenation 
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(3) 

(4) 

(5) 


In  our  previous  report,26  we  followed  an  approach  used  to  measure  regional  cerebral  blood  flow  (rCBF)  with 
diffusible  radiotracers,  as  originally  developed  by  Kety31  in  the  1950’s,  to  develop  our  model.  By  applying 
Fick’s  principle  and  defining  y  as  the  ratio  of  Hb02  concentration  changes  in  the  vascular  bed  to  that  in 
veins,  which  is  equal  to  (AHb02vascubture)/(  AHb02vein),  we  could  get  Eq.  (8)  as  follows: 

AHb02vasculaturc(/)  =  jHJl-expt-yr/Y)]  =  A,[  l-exp(-t/x,)]  (8) 

where  y  is  the  vasculature  coefficient  of  the  tumor.  Ho  is  the  arterial  oxygenation  input  and  /  is  the  blood 
perfusion  rate. 

If  a  tumor  has  two  distinct  perfusion  regions,  and  the  measured  signal  is  from  both  regions  (Figure  2),  then  it 
is  reasonable  that  we  will  get  two  different  blood  perfusion  rates  f\  and  /2,  two  different  vasculature 
coefficients  y,  and  y>,  or  all  four.  Therefore,  Eq.  (8)  can  be  modified  to  double  exponential  expression  as 

AHbO2vasculature(0  =  tfHJl-exp  (-/ir/tf)]  +  ftHJl-expC-/^)] 

=  A,[l-exp(-f/x,)]  +  A2[l-exp(-f/x2)]  (9) 

where  f\  and  y,  are  the  blood  perfusion  rate  and  the  vasculature  coefficient,  respectively,  in  region  1,  /2  and 
yi  are  the  same  for  region  2,  and  A!  =  ^Hq,  A2  =  Xi  =  /,/yi,  x2  =  /2/^.  Then,  if  Ai,  A2,  Xj,  and  x2  are 
determined  from  our  measurements,  we  can  obtain  the  ratios  of  two  vasculature  coefficients  and  the  two 
blood  perfusion  rates: 


Zl  =  A,  /i  _  AJAi 

Yl  ^2  f 2  *{1*2 


(10) 
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With  these  two  ratios,  we  will  be  able  to  understand  more  about  tumor  physiology  such  as  vasculature  and 
blood  perfusion.  In  this  report,  we  used  three-channel  NIR  spectroscopy  to  explore  the  tumor  vascular 
heterogeneity  by  obtaining  At,  A2,  Ti,  and  t2  from  three  different  locations  on  tumors. 

3.  RESULTS 

Figure  3(a)  shows  the  A[Hb02]  obtained  from  three  detectors  in  a  small  breast  tumor  (1.54  cm3),  with  a 
source-detector  separation  of  1  cm  for  detector  #1,  1.3  cm  for  detector  #2,  and  1.8  cm  for  detector  #3, 
respectively.  Detectors  #1  and  #2  are  placed  on  the  side  of  tumor  in  a  reflection  mode  while  detector  #3  was 
located  across  of  light  source  in  a  transmittance  mode.  The  measurement  uncertainties  are  shown  at  only 
discrete  locations.  After  5  min  of  air  breathing  measurement  as  a  baseline,  inhaled  gas  was  switched  from  air 
to  carbogen,  which  caused  a  sharp  increase  in  A[Hb02]  (p  <  0.0001  after  1  min  from  gas  switch)  followed  by 
a  further  slow,  gradual,  but  significant,  increase  over  the  next  15  min  (P  <  0.0001). 

The  rising  part  of  A[Hb02]  from  detector  #1,  #2,  and  #3  after  gas  switch  from  air  to  carbogen  and  the 
fitted  curves  are  shown  in  Figure  3(b)  -  3(d),  respectively.  We  used  single-exponential  (Eq.  8)  and  double¬ 
exponential  expressions  (Eq.  9)  to  fit  the  data  of  rising  portion  of  A[Hb02]  and  it  appears  that  double¬ 
exponential  expression  gives  a  slightly  better  fit,  as  is  confirmed  by  the  respective  R  values  [0.99  versus  0.98 
(not  shown  in  the  table)].  Time  constants  and  amplitudes  from  three  fitted  curves  are  summarized  in  the  Table 
1  to  describe  tumor  vascular  structure  and  blood  perfusion.  By  manipulating  time  constants  and  amplitudes, 
we  can  get  physiologically  meaningful  parameters  such  as  Vi/y2=  Ai/A2  and  /i//2  =  (Ai/A2)/(Ti/t2)  (Eq.  10). 
When  Yi /y2  is  close  to  1,  it  implies  that  the  measured  signal  results  equally  from  both  region  1  and  2,  and  if 
Y1/Y2  <  1,  then  the  measured  signal  results  more  from  region  1  than  region  2  [Figure  3(b),  3(c)  and  3(d)].  All 
the  ratios  of /i/f2  in  Figure  3(b)  -  3(d)  are  larger  than  1,  which  indicate  the  blood  perfusion  rate  in  region  1  is 
much  higher  that  that  in  region  2. 
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Figure  4  depicts  A[Hb02]  obtained  from  another  large  breast  tumor  (16.6  cm3).  The  fast  and  significant 
increases  of  A[Hb02]  after  switching  gas  from  air  to  carbogen  are  similar  to  those  from  tumor  in  Figure  3. 
However,  double-exponential  fitting  shows  much  better  fitness  compared  to  single-exponential  expression  in 
this  tumor  (R=0.79~0.84  from  single  exponential  model  fitting)  and  also  this  tumor  has  several  different 
characteristics  compared  to  tumor  in  Figure  3  as  can  be  seen  from  Table  2.  First,  this  tumor  has  very  different 
values  of  time  constants  between  detectors,  while  the  tumor  in  Figure  3  has  close  values  of  time  constants 
between  detectors.  Fast  time  constants  in  Figure  3  are  0.25,  0.23,  and  0.26  min,  and  slow  time  constants  in 
Figure  3  are  9.16,  8.06,  and  8.52  min,  which  do  not  have  a  significant  differences  between  detectors  (p  >  ), 
but  tumor  in  Figure  4  shows  very  different  values  of  time  constant  between  detectors  (p<0.001  between 
detector  #  and  #  and  p<  between  detector  #  and  #  ).  Fast  time  constants  in  Figure  4  are  0.24, 

0.30,  and  0.27  min,  and  slow  time  constants  in  Figure  4  are  8.3,  9.9,  and  7.0  min.  In  addition  to  time 
constants,  the  Y1/Y2  and  /i//2  are  also  very  different  from  those  in  the  Table  1.  Figures  4(b)  and  4(d)  have  Y1/Y2 
>1  while  all  the  detectors  in  Figure  3  show  <1.  which  implies  that  signals  from  detectors  #1  and  #3  in 
Figure  4  are  mainly  from  region  1  rather  than  region  2.  The  ratios  of  perfusion  rate  in  Figure  4 ,/i//2,  are  also 
much  greater  than  those  from  tumor  in  Figure  3. 

4.  DISCUSSION 

We  have  seen  tumor  heterogeneity  in  terms  of  tissue  p02  using  19F  MR  p02  mapping  via  reports  from  Dr. 
Mason’s  research  group.14'15'16  19F  MR  p02  mapping  technique  clearly  showed  that  tumor  tissue  has  a  wide 
distribution  of  p02,  and  also  different  responses  to  respiratory  challenges  even  within  a  tumor.  We  know  that 
tissue  p02  is  mainly  decided  by  the  balance  between  the  supply  of  oxygen  from  blood  vessels,  and  the  oxygen 
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consumption  by  tissue  cells.  With  assumption  of  constant  oxygen  consumption  during  gas  switch  from  air  to 
carbogen,  tumor  tissue  p02  will  be  mainly  dependent  of  blood  oxygen  level  changes. 

NIRS  is  a  portable,  low  cost,  real  time  measurement  system  that  can  measure  changes  of  blood  oxygen  level 
by  using  two  wavelengths.  We  have  used  a  NIRS  system  that  has  one  light  source  and  one  detector  so  that  we 
could  have  global  measurement  of  A[Hb02]  and  A[Hbtotai]  in  tumor  during  respiratory  challenges.  Through 
our  previous  experiments,  we  found  that  most  of  tumors  have  bi-phasic  behavior  of  A[Hb02]  increase  (rapid 
increase  followed  by  a  slow  and  gradual  increase)  after  switching  the  gas  from  air  to  carbogen  or  air  to  100% 
oxygen.  To  explain  this  bi-phasic  behavior,  we  formed  a  hypothesis  that  the  bi-phasic  behavior  of  A[Hb02] 
increase  during  carbogen  or  oxygen  inhalation  results  from  two  different  regions  in  tumors  in  terms  of  blood 
perfusion  rate  and  vascular  structure,  and  we  developed  a  mathematical  model  to  get  physiologically 
meaningful  parameters  from  experimental  data. 

The  region  in  tumor  that  has  fast  blood  perfusion  rate,/i,  and  high  ratio  of  A[Hb02]vasculature  and  A[Hb02]vem, 
Yi,  is  called  region  1,  while  region  2  represents  the  region  that  has  slow  blood  perfusion  rate,/2,  and  low  ratio 
of  A[Hb02]vascuIature  and  A[Hb02]vein,  y2-  Thus,  the  ratio  of  Yi  and  Y2  can  be  simply  thought  of  as  a  ratio  of 
region  1  volume  to  region  2  volume.  In  other  words,  when  a  tumor  has  Y1/Y2  «  1,  it  implies  that  the  tumor 
consists  of  50%  well-perfused  region  and  50%  poorly  perfused  region.  If  a  tumor  has  Y1/Y2  <1  or  Y1/Y2  >1*  this 
means  the  measured  signal  results  mostly  from  poorly  perfused  region,  or  well-perfused  region,  respectively. 
Therefore,  by  looking  at  these  two  ratios  from  Table  1  and  2,  we  can  find  the  differences  between  tumor  in 
Figure  3  and  tumor  in  Figure  4  in  terms  of  blood  perfusion  rate  and  proportion  of  either  well-perfused  region 
or  poorly  perfused  region.  For  instance,  a  large  proportion  of  tumor  in  Fig.  3  consists  of  well-perfused  region 
(Y1/V2  =  0-96  ~  2.23)  while  tumor  in  Figure  4  is  mostly  composed  of  poorly  perfused  region  (Y1/Y2  =  0.135  ~ 
0.343).  In  the  case  of  tumor  in  Figure  3,  the  blood  perfusion  rate  in  the  well-perfused  region  is  much  faster 
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than  in  the  poorly  perfused  region,  while  the  ratio  of  blood  perfusion  rate  between  well-perfused  region  and 
poorly  perfused  region  from  tumor  in  Figure  4  (fx/f2  =  5.13-  18.33)  is  not  as  high  as  tumor  in  Figure  3  (f\/f2  = 
3 1.58 -63.75). 

5.  CONCLUSION 

We  applied  a  three  channel  NIRS  system  to  show  the  heterogeneity  of  tumor  vasculature  in  terms  of  blood 
perfusion  rate  and  vascular  coefficients.  All  detectors  on  tumor  showed  a  sharp  increase  of  A[Hb02]  that  was 
followed  by  slow  and  gradual  increase.  However,  the  differences  between  the  signals  from  three  detectors 
could  be  found  by  applying  our  mathematical  model  to  fit  the  A[Hb02]  data  during  carbogen  inhalation. 
Although  all  signals  from  three  detectors  in  both  Figures  3  and  4  showed  bi-phasic  behavior  of  A[Hb02],  their 
time  constants  and  ratio  of  Yi/y2  and  f\lf2  revealed  that  tumor  vasculature  is  not  homogeneous  and,  indeed,  it  is 
heterogeneous  with  respect  to  blood  perfusion  rate  and  proportion  of  well-perfused  region  in  tumor.  Our 
future  study  will  be  conducted  to  prove  our  hypothesis  on  the  bi-phasic  behavior  of  A[Hb02]  during  carbogen 
inhalation  and  to  learn  more  about  meaning  of  time  constant,  amplitude,  Yi/y2,  and /t//2  from  our  mathematical 
model. 
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FIGURE  CAPTION 

Figure  1.  Schematic  experimental  setup  of  animal  experiment.  One  of  three  detectors  is  placed  opposite  of  the 
light  source  to  detect  light  in  transmission  mode,  and  the  other  two  were  located  in  the  semi-reflection  mode. 

Figure  2.  Schematic  diagram  of  two-region  tumor  model.  If  tumor  has  two  distinct  regions  in  terms  of 
perfusion  rate,  the  detector  will  get  signal  from  both  regions,  and  the  signal  can  be  fitted  with  double 
exponential  expression. 

Figure  3.  Dynamic  changes  of  [Hb02]  from  three  detectors  in  a  rat  breast  (volume:  1.5  cm3)  tumor  A.  Dotted 
vertical  lines  mark  the  point  when  the  gas  was  changed  [Fig.  3(a)].  The  rising  part  of  A[Hb02]  from  detectors 
was  fitted  using  both  single-exponential  and  double-exponential  expressions  and  Fig.  3(b)-3(d)  are  from 
detector  #l-#3,  respectively. 

Figure  4.  Dynamic  changes  of  [Hb02]  from  three  detectors  in  a  rat  breast  (volume:  16.6  cm3)  tumor  B.  Dotted 
vertical  lines  mark  the  point  when  the  gas  was  changed  [Fig.  4(a)].  Figures  4(b)-4(d)  show  the  rising  part  of 
A[Hb02]  after  gas  switch  from  air  to  carbogen,  and  fitted  curves  in  detector  #l-#3,  respectively. 

Figure  5.  Dynamic  changes  of  [Hb02]  from  three  detectors  in  a  rat  breast  (volume:  20.6  cm3)  tumor  C.  Dotted 
vertical  lines  mark  the  point  when  the  gas  was  changed  [Fig.  5(a)].  Figures  5(b)-5(d)  show  the  rising  part  of 
A[Hb02]  after  gas  switch  from  air  to  carbogen,  and  fitted  curves  in  detector  #l-#3,  respectively. 
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Table  1.  Summary  of  vascular  oxygen  dynamics  determined  at  the  three  detectors  from  Figure  3  for  Tumor  A 


Double-Exponential  fitting  AHb02  -  A^l-expC-t/xO]  +  A2[l-exp(-t/x2)] 

Parameters 

Detector  #1  (Fig.  3b) 

Detector  #2  (Fig.  3c) 

Detector  #3  (Fig.  3d) 

Separation:  d  (cm) 

1.0 

1.3 

1.8 

A,  (mM/DPF) 

0.0069  ±  0.0003 

0.0055  1  0.0002 

0.004  1  0.0002 

Xi  (min) 

0.25  1  0.04 

0.23  1 0.04 

0.26 1 0.05 

A2  (mM/DPF) 

0.041  ±0.001 

0.032  1  0.0003 

0.024  1  0.0003 

x2  (min) 

9.16±0.37 

8.06 1 0.27 

8.52  1 0.36 

X2 

0.00047 

0.0003 

0.00021 

R 

0.99 

0.99 

0.99 

Yi/Y2=  A,/A2 

0.17  ±0.01 

0.17  ±0.01 

0.1710.01 

X,/X2 

0.027  ±  0.005 

0.029  1  0.006 

0.031  ±0.007 

/i//2  =  (A,/A2)/(t,/t2) 

6.3  ±  1.9 

5.912 

5.512 

69 


Table  2.  Summary  of  vascular  oxygen  dynamics  determined  at  the  three  detectors  from  Figure  4  for  Tumor  B 


Double-Exponential  fitting  AHb02  =  Ajfl-expO-t/Xi)]  +  A2[l-exp(-t/T2)] 

Parameters 

Detector  #1  (Fig.  4b) 

Detector  #2  (Fig.  4c) 

Detector  #3  (Fig.  4d) 

Separation:  d  (cm) 

1.5 

2.5 

2.8 

A,  (mM/DPF) 

0.0370  ±  0.0004 

0.0125  ±  0.0002 

0.0134  ±  0.0002 

X!  (min) 

0.24  ±0.01 

0.30  ±  0.02 

0.27  ±0.01 

A2  (mM/DPF) 

0.020  ±0.001 

0.013  ±  0.0004 

0.006  ±  0.0001 

x2  (min) 

8.3  ±  0.7 

9.9  ±0.8 

7.0  ±0.6 

x2 

0.0045 

0.001 

0.00046 

R 

0.95  0.79 

0.96  0.82 

0.95  0.84 

Yi/Y2=  A,/A2 

1.85  ±0.07 

0.96  ±  0.05 

2.23  ±  0.08 

Xf  fh 

0.029  ±  0.004 

0.03  ±  0.004 

0.039  ±0.005 

/,//2  =  (A1/A2)/(t1/x2) 

64  ±9 

32  ±5 

56  ±9 
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Table  3.  Summary  of  vascular  oxygen  dynamics  taken  at  the  three  detectors  from  Figure  5  for  Tumor  C 


Double-Exponential  fitting  AHb02  =  Ai[l-exp(-t/ii)]  +  A2[l-exp(-t/x2)] 

Parameters 

Detector  #1  (Fig.  4b) 

Detector  #2  (Fig.  4c) 

Detector  #3  (Fig.  4d) 

Separation:  d  (cm) 

2 

2.5 

3 

A,  (mM/DPF) 

0.01  ±0.0003 

0.0072  ±  0.0001 

0.0023  ±  0.0001 

Tj  (min) 

0.27  ±  0.03 

0.43  ±  0.02 

0.16  ±0.04 

A 2  (mM/DPF) 

0.033  ±  0.0003 

0.021  ±0.002 

0.017  ±  0.0001 

t2  (min) 

4.7  ±0.1 

23  ±3 

6.1  ±0.2 

x2 

0.00093 

0.00019 

0.00031 

R 

0.99  0.97 

0.99  0.91 

0.99  0.98 

P* 

< 

< 

II 

4? 

0.30  ±0.01 

0.34  ±  0.03 

0.1 35  ±0.009 

tl/X2 

0.058  ±  0.007 

0.019  ±0.003 

0.026  ±  0.007 

/,//2  =  (A,/A2)/(t,/t2) 

5.2  ±0.7 

18  ±  4 

5  ±  1 
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Appendix  3 


Dynamic  response  of  breast  tumor  oxygenation  to 
hyperoxic  respiratory  challenge  monitored  with 
three  oxygen-sensitive  parameters 


Yueqing  Gu,  Vincent  A.  Bourke,  Jae  G.  Kim,  Anca  Constantinescu,  Ralph  P.  Mason, 
and  Hanli  Liu 


The  simultaneous  measurement  of  three  oxygen-sensitive  parameters  [arterial  hemoglobin  oxygen  sat¬ 
uration  (Sa02),  tumor  vascular-oxygenated  hemoglobin  concentration  ([HbO  J),  and  tumor  oxygen  ten¬ 
sion  (p02)]  in  response  to  hyperoxic  respiratory  challenge  is  demonstrated  in  rat  breast  tumors.  The 
effects  of  two  hyperoxic  gases  [oxygen  and  carbogen  (5%  C02  and  95%  02)]  were  compared,  by  use  of  two 
groups  of  Fisher  rats  with  subcutaneous  13762NF  breast  tumors  implanted  in  pedicles  on  the  foreback. 
Two  different  gas-inhalation  sequences  were  compared,  i.e.,  air- carbogen-air- oxygen-air  and  air- 
oxygen-air-  carbogen-air .  The  results  demonstrate  that  both  of  the  inhaled,  hyperoxic  gases  signifi¬ 
cantly  improved  the  tumor  oxygen  status.  All  three  parameters  displayed  similar  dynamic  response  to 
hyperoxic  gas  interventions,  but  with  different  response  times:  the  fastest  for  arterial  Sa02,  followed  by 
biphasic  changes  in  tumor  vascular  [Hb02],  and  then  delayed  responses  for  p02.  Both  of  the  gases 
induced  similar  changes  in  vascular  oxygenation  and  regional  tissue  p02  in  the  rat  tumors,  and  changes 
in  [HbOJ  and  mean  p02  showed  a  linear  correlation  with  large  standard  deviations,  which  presumably 
results  from  global  versus  local  measurements.  Indeed,  the  p02  data  revealed  heterogeneous  regional 
response  to  hyperoxic  interventions.  Although  preliminary  near- infrared  measurements  had  been  dem¬ 
onstrated  previously  in  this  model,  the  addition  of  the  p02  optical  fiber  probes  provides  a  link  between 
the  noninvasive  relative  measurements  of  vascular  phenomena  based  on  endogenous  reporter  molecules, 
with  the  quantitative,  albeit,  invasive  p02  determinations.  ©  2003  Optical  Society  of  America 
OCIS  codes:  170.1470,  170.3660,  170.4580,  120.3890,  120.1880,  230.2090. 


1.  Introduction 

It  is  widely  recognized  that  hypoxic  regions  in  solid 
tumors  may  limit  the  efficacy  of  nonsurgical  therapy, 
including  radiotherapy,  photodynamic  therapy,  and 
chemotherapy.1-4  Many  adjuvant  interventions 
have  been  tested,  including  simple  strategies  such  as 
breathing  hyperoxic  gases.5"7  However,  a  meta¬ 
analysis  of  some  10,000  patients  showed  only  a  mod¬ 
est  benefit,  and  this  benefit  was  restricted  to  specific 
tumor  types.8  It  is  thought  that  the  failure  of  such 
interventions  was  largely  due  to  the  inability  to  iden- 
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tify  those  patients  who  would  benefit.  Indeed,  there 
is  growing  emphasis  on  tailoring  therapy  to  the  indi¬ 
vidual  characteristics  of  each  patient’s  tumor.  Fur¬ 
thermore,  carbogen  (5%  C02  and  95%  02)  and  oxygen 
have  been  used  on  experimental  tumors  in  animals  as 
well  as  on  clinical  trials  in  patients  for  many 
years.910  But  the  therapeutic  benefits  of  the  two 
kinds  of  respiratory  hyperoxic  gases  are  diverse,  de¬ 
pending  on  the  tumor  types  and  individuals.11"13 
Accordingly,  accurate  assessment  of  tumor  oxygen¬ 
ation  at  various  stages  of  tumor  growth  and  in  re¬ 
sponse  to  interventions  may  provide  a  better 
understanding  of  tumor  development  and  may  serve 
as  a  prognostic  indicator  for  treatment  outcome,  po¬ 
tentially  allowing  therapy  to  be  tailored  to  individual 
characteristics. 

Various  techniques  have  been  developed  to  mea¬ 
sure  oxygen  tension  (p02)  or  vascular  oxygenation  of 
tumors.14  Many  methods  are  invasive,  and  those 
requiring  biopsy  preclude  dynamic  investigations. 
Optical  techniques  based  on  light  absorption  of  en¬ 
dogenous  chromophores,  e.g.,  near-infrared  (NIR) 
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spectroscopy  of  oxygenated  and  deoxygenated  hemo¬ 
globin,  are  entirely  noninvasive  and  allow  real-time 
monitoring  of  tumor  vascular  oxygenation.15-17 
However,  NIR  has  limited  spatial  resolution,  and  it 
remains  to  be  determined  whether  vascular  oxygen¬ 
ation  is  related  to  therapeutic  outcome.  Hitherto, 
quantitative  p02  has  been  shown  to  have  prognostic 
value,18-21  but  p02  represents  a  balance  between  ox¬ 
ygen  delivery  and  consumption.  Thus,  we  seek  to 
explore  the  interplay  of  vascular  and  tissue  oxygen¬ 
ation.  Electrodes  have  been  used  widely  to  study 
tumor  oxygen  dynamics  with  respect  to 
interventions,22-24  but  they  are  generally  limited  to  a 
single  location  and  small  probes  can  be  fragile.  We 
have  ourselves  recently  shown  a  correlation  between 
p02  and  AHb02  in  some  tumors,  but  we  noted  dis¬ 
tinct  heterogeneity,  and  thus,  the  global  NIR  mea¬ 
surements  were  not  always  related  to  local  p02.25 
Multiple  fiber-optic  probes  may  be  inserted  into  a 
tumor,26-28  and  we  have  now  investigated  correlation 
between  NIR  measurements  and  multiple  (three)  si¬ 
multaneous  p02  measurements. 

We  now  report  simultaneous  measurements  of 
three  oxygen-related  parameters,  i.e.,  arterial  hemo¬ 
globin  oxygen  saturation,  Sa02;  tumor  oxygenated 
hemoglobin  concentration,  [Hb02];  and  tumor  oxygen 
tension,  p02,  to  assess  dynamic  responses  of  rat 
breast  tumors  to  hyperoxic  gases .  Changes  in  tumor 
vascular  [Hb02]  were  measured  by  NIR  spectroscopy 
(NIRS)  using  a  photon-migration,  frequency-domain 
device;  changes  in  regional  p02  were  monitored  by  a 
fluorescence-quenched,  oxygen-sensing,  fiber-optic 
system  (FOXY);  the  arterial  Sa02  values  were  re¬ 
corded  with  a  fiber-based,  pulse  oximeter. 

2.  Materials  and  Methods 

A.  Near-Infrared  Spectroscopy  System  for  Measurement 
of  Changes  in  [Hb02] 

NIR  light  (700  to  900  nm)  has  considerable  tissue 
penetration  depth  (several  centimeters)  and  permits 
in  vivo  sampling  of  large  tissue  volumes  (e.g.,  human 
breast,  brain,  skeletal  muscle,  or  tumors),  since  pho¬ 
ton  transport  in  tissue  is  dominated  by  scattering 
rather  than  by  absorption.  Absorption  of  NIR  light 
by  the  oxygenated  and  the  deoxygenated  hemoglobin 
chromophores  may  be  used  to  determine  hemoglobin 
oxygenation  and  blood  concentration  changes.  As 
described  in  detail  previously,1625  a  homodyne 
frequency-domain  system  (NIM,  Philadelphia,  Penn¬ 
sylvania)  was  used  to  monitor  the  global  changes  in 
oxygenated  and  deoxygenated  hemoglobin  concentra¬ 
tions,  A[Hb02]  and  A[Hb],  respectively,  in  rat  breast 
tumors  in  response  to  variations  in  inhaled  gas. 
Briefly,  the  light  from  two  NIR  laser  diodes  (758  nm 
and  785  nm)  was  coupled  into  a  bifurcated  fiber  bun¬ 
dle  and  illuminated  on  the  tumor,  and  the  transmit¬ 
ted  light  was  collected  and  propagated  to  a 
photomultiplier  tube  (Fig.  1).  The  fiber  bundles 
were  placed  on  the  surface  of  the  tumors  in  a  trans¬ 
mittance  mode  parallel  to  the  body  of  the  rat.  The 
fiber  tips  touched  firmly  on  the  skin  (without  com- 


FOXY  System 


Fig.  1.  Experimental  setup  for  simultaneous  oximetry.  The 
3-mm-diameter  fiber  bundles  of  the  NIRS  system  deliver  and  de¬ 
tect  the  laser  light  through  the  tumor  in  transmittance  geometry. 
PMT  represents  a  photomultiplier  tube.  I/Q  is  an  in-phase  and 
quadrature  phase  demodulator  for  retrieving  amplitude  and  phase 
information.  The  FOXY  system  comprises  three  fiber-optic 
oxygen-sensing  probes  that  are  inserted  into  different  regions  of 
the  tumor.  The  pulse  oximeter  probe  is  placed  on  the  hind  foot  of 
the  rat. 


pression)  in  the  middle  parts  of  the  tumors,  providing 
optimal  geometry  to  interrogate  deep  tumor  tissue. 

Based  on  modified  Beer-Lambert’s  law,29  changes 
in  oxygenated  and  deoxygenated  hemoglobin  concen¬ 
trations,  A[Hb02]  and  A[Hb],  due  to  respiratory  in¬ 
tervention  were  derived  from  the  measured 
amplitudes  at  the  two  wavelengths  and  calculated 
with  the  following  equations25: 


A[Hb02]  = 


10.63  log(^ 


758 


+  14.97  log)^5 


785 


(l) 


8.95  log) 


A[Hb]  = 


\758  I A  \786 


(2) 


where  AB  and  AT  are  the  baseline  and  transient  am¬ 
plitudes  measured  from  the  NIR  system,  respective¬ 
ly;  d  is  the  source- detector  separation;  the  unit  for 
both  A[Hb02]  and  A[Hb]  is  millimolar  per  differential 
path-length  factor  (DPF);  and  the  DPF  is  for  tumor 
tissues.  As  demonstrated  in  our  previous  study, 
normalization  of  A[Hb02]  and  A[Hb]  to  their  maximal 
values  can  eliminate  the  effects  of  d  and  DPF  on  the 
results.25 


B.  Fiber-Optic  Oxygen-Sensing  System  for  Measurement 
of  Changes  in  p02 

Regional  p02  in  tumors  was  monitored  with  a  mul¬ 
tichannel,  fiber-optic,  oxygen-sensing  system  (FOXY, 
Ocean  Optics,  Inc.,  Dunedin,  Florida).30  Three 
fluorescence-quenched,  optical  fiber  probes  (AL300, 
tip  diameter  410  |xm)  were  inserted  into  different 
regions  of  the  tumors  (Fig.  1).  Probes  were  posi¬ 
tioned  so  that  at  least  one  was  in  a  poorly  oxygenated 
region  (low  baseline  p02)  and  at  least  one  in  a  well- 
oxygenated  region  (high  baseline  p02).  If  necessary, 
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the  probes  were  gently  moved  through  the  tumor 
until  such  diverse  regions  were  located.  In  some 
cases,  the  mean  p02  derived  from  the  three  individ¬ 
ual  measurements  is  presented.  Although  this  is  a 
commercial  system,  few  details  have  been  published 
previously,31  and  no  applications  to  in  vivo  tumor 
oximetry  have  been  published  to  our  knowledge. 
Light  from  a  pulsed  blue  LED  (475  nm)  was  coupled 
into  one  branch  of  a  bifurcated  optical  fiber  bundle 
and  propagated  to  the  probe  tip.  The  distal  end  of 
the  probe  is  coated  with  a  thin  layer  of  a  hydrophobic 
solgel  material,  where  an  oxygen-sensing  ruthenium 
complex  is  effectively  trapped.  Illumination  of  the 
ruthenium  complex  causes  fluorescence  at  —600  nm. 
If  the  excited  ruthenium  complex  encounters  an  ox¬ 
ygen  molecule,  the  excess  energy  is  transferred  to  the 
oxygen  molecule  in  a  nonradiative  transition,  de¬ 
creasing  or  quenching  the  fluorescence  signal.  The 
degree  of  quenching  correlates  with  the  oxygen  con¬ 
centration,  and  hence,  p02. 

The  fluorescence  response  of  the  ruthenium  crystal 
complex  is  highly  temperature  dependent,  so  to  ac¬ 
complish  probe  calibration  it  was  necessary  to  stream 
gases  of  known  oxygen  concentrations  (100%,  20.9%, 
10%,  2%,  and  0%)  through  a  cylindrical  water  jacket 
heated  to  37  °C.  Calibration  curves  were  automati¬ 
cally  calculated  by  means  of  the  vendor-supplied  soft¬ 
ware,  with  the  second-order,  polynomial  calibration: 

y  =  1  +  #i[0]  +  K2[Of  (3) 

where,  I0  is  the  fluorescence  intensity  at  zero  oxygen 
concentration  (nitrogen),  /  is  the  measured  intensity 
of  fluorescence  at  a  pressure  of  oxygen,  [O]  represents 
the  oxygen  concentration  (related  to  p02),  K1  and  K2 
are  the  first-  and  the  second-order  coefficients  and  are 
automatically  supplied  by  the  curve-fitting  routine 
from  the  calibration  measurements. 

C.  Pulse  Oximeter  for  Measurement  of  Arterial  Sa02 

Arterial  Sa02  of  the  breast-tumor-bearing  rats  was 
also  monitored  with  a  fiber-optic  pulse  oximeter 
(Nonin  Medical,  Inc.,  Plymouth,  Minnesota)  placed 
on  the  hind  foot  of  the  rats.  The  system  consisted  of 
two  optical  fibers  used  for  delivering  and  receiving 
the  light.  The  tips  were  placed  on  either  side  of  the 
foot  in  transmission  mode. 

D.  Animal  Model 

Mammary  adenocarcinomas  13762NF  (originally  ob¬ 
tained  from  the  Division  of  Cancer  Therapeutics, 
NIH,  Bethesda,  Maryland)  were  implanted  in  skin 
pedicles32  on  the  foreback  of  adult  female  Fisher  344 
rats  (—150  g).  Once  the  tumors  reached  1-2  cm  di¬ 
ameter,  rats  were  anesthetized  with  150- jxl  ketamine 
hydrochloride  (100  mg/ml,  i.p.)  and  maintained  un¬ 
der  general  gaseous  anesthesia  with  1.3%  isoflurane 
in  air  (1  dm3/min).  Body  temperature  was  main¬ 
tained  at  37  °C  by  a  warm  water  blanket.  Tumors 
were  shaved  to  improve  optical  contact  for  transmit¬ 
ting  NIR  light.  The  tumor  diameters  along  the 


air  carbogen  air 


Fig.  2.  Time  profile  of  the  three  oxygen-sensitive  parameters,  i.e., 
the  normalized  changes  of  tumor  AfHbOJ,  the  mean  changes  of 
tumor  Ap02,  and  the  arterial  Sa02  with  respect  to  carbogen 
breathing  in  a  representative  13762NF  rat  breast  tumor  (No.  1, 3.2 
cm3). 


three  major  orthogonal  axes  (a,  6,  c)  were  measured 
with  calipers  and  volume  estimated  with  an  ellipsoid 
approximation  with  the  formula:  V  =  (n/6)(abc). 

Two  groups  of  rats  (n  =  7  in  each  group)  were  used 
to  compare  the  effects  of  carbogen  and  oxygen  on 
vascular  oxygenation  of  breast  tumors.  Group  1  ex¬ 
perienced  the  gas-inhalation  sequence  of  air- 
carbogen-air- oxygen-air.  Group  2  was  exposed  to 
the  reverse  sequence  of  air-  oxygen-air-  carbogen- 
air.  Each  gas  was  maintained  for  20  min.  In  addi¬ 
tion,  the  FOXY  p02  probes  were  applied  to  five  rats 
from  Group  1,  and  the  dynamics  of  the  three  oxygen- 
related  parameters  were  measured  simultaneously. 

3.  Results 

A.  Dynamic  Responses  of  Three  Oxygen-Related 
Parameters  to  Carbogen  Intervention 

Typical  time  profiles  of  the  normalized  A[Hb02], 
mean  Ap02,  and  Sa02  in  response  to  carbogen  inter¬ 
vention  are  shown  for  a  representative  13762NF 
breast  tumor  (No.  1,  3.2  cm3)  in  Fig.  2.  When  the 
inspired  gas  was  switched  from  air  to  carbogen,  the 
Sa02  readings  increased  rapidly  and  significantly 
from  the  baseline  value  of  85%  to  the  maximum  of 
100%  within  2.5  minutes  (p  <  0.0001).  The  normal¬ 
ized  A[Hb02]  showed  a  sharp  initial  rise  in  the  first 
minute  (p  <  0.0001)  followed  by  a  slower,  gradual, 
but  further  significant  increase  over  the  next  19  min 
(p  <  0.001).  Mean  Ap02  increased  rapidly  by  ap¬ 
proximately  50  Torr  within  8  min  (p  <  0.0005)  and 
also  continued  a  slower  and  gradual  increase  over  the 
next  12  min  (p  <  0.005).  Return  to  breathing  air 
produced  a  significant  decline  for  all  three  signals 

(p  <  0.0001). 

Sa02  and  p02  displayed  a  single-phase  dynamic 
behavior  in  response  to  carbogen  intervention, 
whereas  A[Hb02]  showed  an  apparent  biphasic  re¬ 
sponse.  These  dynamics  may  be  characterized  by 
time  constants  of  a  single-exponential  response.  For 
the  tumor  in  Fig.  3,  Sa02  had  the  fastest  response, 
with  a  time  constant  of  T(Sa02)  =  1.1  ±  0.2  min  ( R  = 
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Tim*  (min) 

Fig.  3.  Dynamic  responses  of  the  three  oxygen-sensitive  param¬ 
eters  to  carbogen  intervention  in  a  rat  breast  tumor  (No.  1,  3.2 
cm3).  Single-exponential  curve  fitting  yielded  Sa02  = 
0.204(1  -  exp[-(*  -  20.02)/l.l]}  +  0.85  (. R  =  0.93),  A[Hb02]  = 
0.655(1  -  exp[-(*  -  20.36)/2.59]}  +  0.125  (R  =  0.89),  and 
Ap02  =  42.68(1  -  exp[-(f  -  21.01)/4.56]}  +  16.66  (. R  =  0.98); 
biexponential  fitting  resulted  in  A[Hb02]  =  0.373(1  -  exp[-(f  - 
20.36)/0.61]}  +  0.648(1  -  exp[-(f  -  20.36)/21]}  (R  =  0.97). 


0.93),  followed  by  [Hb02]  with  T(A[Hb02])  =  2.59  ± 
0.06  min  (R  =  0.89),  whereas  Ap02  yielded  the  slow¬ 
est  response  T(Ap02)  =  4.56  ±  0.06  min  (R  =  0.98). 
Time  constants  for  Group  1  are  listed  in  Table  1.  In 
every  case  T(Sa02)  <  T(A[Hb02])  <  T(Ap02),  based  on 
the  single-exponential  fitting.  No  apparent  relation 
between  the  time  constant  and  the  tumor  volume  was 
observed. 

It  is  clear  that  the  response  of  AHb02  is  not  well 
represented  by  a  single  exponential,  and  thus,  a 
double-exponential  expression  with  two  time  con¬ 
stants,  Tj  and  t2,  was  also  used  (Fig.  3).  Comparison 
between  the  biexponential  fitting  for  A[Hb02]  and  the 
single-exponential  results  for  both  Sa02  and  Ap02  in 
the  first  five  rat  tumors  (Table  1)  shows  that  the  time 
constants  of  Sa02  (—1.2  ±  0.4  min)  are  similar  to 
those  of  the  first  phase  of  A[Hb02]  (—0.5  ±  0.2  min), 
whereas  the  second  phase  is  longer  and  highly  vari¬ 


able  (—14  ±  11  min).  No  significant  correlations 
were  found  between  any  of  the  time  constants  in 
Table  1. 

Time  delay,  td ,  between  the  time  when  the  gas 
intervention  was  initiated  and  the  time  when  the 
changes  in  signals  were  detected,  reveals  another 
difference  among  the  three  oxygen-sensitive  param¬ 
eters.  For  tumor  1  (Fig.  2),  the  Sa02  signal  was  the 
first  to  respond  to  the  intervention.  Change  in 
A[Hb02]  was  observed  30  s  later  with  td  -  30  s, 
followed  by  changes  in  Ap02  another  30  s  later  (td  = 
60  s).  Similarly,  when  the  gas  was  returned  from 
carbogen  to  air,  the  Sa02  signal  decreased  immedi¬ 
ately,  followed  by  declines  in  A[Hb02]  and  in  Ap02 
with  td  of  30  and  120  s  later,  respectively.  As  ex¬ 
pected,  changes  in  Sa02  always  preceded  AHb02,  and 
Ap02  occurred  last  for  all  tumors. 

B.  Comparison  of  the  Effects  of  Carbogen  and  Oxygen 
Intervention  on  Tumor  Oxygenation 

Switching  from  air  breathing  to  carbogen  or  oxygen 
produced  similar  changes  in  AHb02  [Fig.  4(a)]. 
However,  the  time  course  was  substantially  different, 
requiring  a  biphasic  exponential  fit  for  carbogen,  but 
a  single  exponential  for  oxygen  [Fig.  4(b)].  For  the 
seven  tumors  in  Group  1,  there  was  no  significant 
difference  (p  >  0.3)  in  the  maximum  magnitude  of 
A[Hb02]  caused  by  carbogen  or  oxygen  interventions 
[Fig.  4(c)]. 

To  examine  the  possible  effect  of  preconditioning 
required  that  Group  2  experience  a  reversed  gas  in¬ 
tervention,  with  exposure  to  oxygen  prior  to  carbogen 
[Fig.  5(a)].  In  this  case,  the  time  constants  of  the 
normalized  tumor  vascular  A[Hb02]  were  now  simi¬ 
lar  for  both  gas  challenges:  indeed,  for  six  of  seven 
tumors,  carbogen  no  longer  induced  the  biphasic  be¬ 
havior.  Figure  5(b)  shows  that  changes  in 
(AHb02)max  were  similar  to  those  in  Group  1,  and 
again,  the  two  gases  did  not  produce  significantly 
different  response  in  (AHb02)max.  This  is  empha¬ 
sized  for  both  Groups  1  and  2  by  a  strong  linear 
correlation  (slope  =  1.16)  between  the  A[Hb02]max 
values  observed  in  response  to  each  of  the  two  con- 


Table  1.  Time  Constants  of  Sa02,  A[Hb02],  and  ApOz  Response  to  Carbogen  and  Oxygen  Intervention  in  the  Breast  Tumors'1 


Breast  Tumors 
Volume  (cm3) 

Single-Exponential  Fitting  of  Sa02,  A[Hb02]  and  ApO 
(Carbogen  Intervention) 

*2 

Double-Exponential  Fitting  for 
A[Hb02]  (Carbogen 
Intervention) 

Single-Exponential 
Fitting  of  A[HbOJ 
(02  Intervention) 

Sa02 

A[Hb02] 

Ap02 

t  (min) 

R 

t  (min) 

R 

t  (min) 

R 

t1  (min) 

t2  (min) 

R 

t  (min) 

R 

No.  1  (3.2) 

1.1  ±  0.2 

0.93 

2.59  ±  0.06 

0.89 

4.56  ±  0.04 

0.98 

0.61  ±  0.03 

21  ±3 

0.97 

0.35  ±  0.01 

0.92 

No.  2  (3.0) 

1.6  ±  0.2 

0.98 

3.40  ±  0.07 

0.91 

4.6  ±  0.1 

0.82 

0.62  ±  0.06 

11  ±  1 

0.96 

0.51  ±  0.01 

0.91 

No.  3  (4.6) 

1.2  ±  0.2 

0.97 

2.12  ±  0.06 

0.76 

2.26  ±  0.02 

0.98 

0.6  ±  0.1 

37  ±3 

0.96 

1.52  ±  0.02 

0.89 

No.  4  (2.6) 

1.9  ±  0.3 

0.98 

2.68  ±  0.05 

0.93 

3.5  ±  0.1 

0.86 

0.12  ±  0.02 

5.2  ±  0.1 

0.98 

1.71  ±  0.03 

0.94 

No.  5  (5.6) 

0.8  ±  0.2 

0.91 

2.68  ±  0.05 

0.74 

4.51  ±  0.02 

0.99 

0.17  ±  0.03 

12.5  ±  0.6 

0.99 

5.49  ±  0.03 

0.98 

No.  6  (1.9) 

0.9  ±  0.2 

0.81 

1.62  ±  0.01 

0.95 

nd 

/ 

0.63  ±  0.08 

2.3  ±  0.1 

0.96 

5.16  ±  0.06 

0.93 

No.  7  (0.72) 

1.0  ±  0.5 

0.95 

3.60  ±  0.03 

0.93 

nd 

/ 

0.61  ±  0.02 

10.5  ±  0.3 

0.98 

3.54  ±  0.03 

0.95 

Mean 

1.2  ±  0.4 

2.7  ±  0.6 

4  ±  1 

0.5  ±  0.2 

14  ±  11 

2.5  ±  2 

“Under  the  inhalation  sequence  of  air-  carbogen-air-  oxygen-air. 
6nd,  not  determined. 
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Fig.  4.  (a)  Time  course  of  changes  in  tumor  vascular  A[HbOJ  for 
a  representative  13762NF  breast  tumor  from  Group  1  (No.  2,  3.0 
cm3)  with  respect  to  altering  inhaled  gas.  (b)  Respective  curve  fits 
for  the  carbogen  and  oxygen  interventions,  (c)  Average  maximum 
values  of  normalized  A[Hb02]  for  the  seven  breast  tumors  in  Group 

L 


secutive  interventions  [Fig.  5(c)].  In  this  case,  non- 
normalized  data  are  shown  for  specific  comparison  of 
the  absolute  A[Hb02]max  produced  by  oxygen  and  car¬ 
bogen  for  each  of  the  tumors. 

C.  Tumor  p02 

The  FOXY  p02  probes  generally  indicated  distinct 
heterogeneity  in  p02.  Moreover,  response  to  the  hy- 
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Fig.  5.  (a)  Dynamic  changes  in  tumor  vascular  A[Hb02]  for  a 
representative  13762NF  breast  tumor  from  Group  2  (No.  9,  2.6 
cm3)  with  gas-inhalation  sequence  reversed  compared  with  Group 

1.  (b)  Average  maximum  values  of  normalized  A[Hb02]  for  Group 

2.  Gas-inhalation  sequence  reversed  compared  with  Group  1. 

(c)  Correlation  between  maximum  A[Hb02]  achieved  with  carbo¬ 
gen  inhalation  versus  that  with  oxygen  {R  =  0.97):  ♦,  carbogen 

first;  A,  oxygen  first. 


peroxic  gas  was  diverse:  those  probes  that  indicated 
apparently  well-oxygenated  regions  usually  showed  a 
large  and  rapid  response,  whereas  those  with  lower 
baseline  p02  often  showed  little  change  [Fig.  6(a)]. 
The  p02  responses  to  the  two  interventions  showed  a 
highly  consistent  behavior  at  each  individual  location 
[Fig.  6(b)].  There  was  also  a  distinct  correlation  be¬ 
tween  the  global  NIR  measurements  and  the  mean 
Ap02  (Fig.  7).  Because  of  heterogeneity  in  regional 
p02,  the  standard  deviations  of  the  mean  p02  values 
were  large. 


air  carbogen  air  oxygen  air 


Fig.  6.  (a)  Time  profiles  of  tumor  Ap02,  measured  with  the  three 
channels  of  the  FOXY  fiber-optic,  oxygen-sensing  system  with  re¬ 
spect  to  different  gas  inhalations  for  breast  tumor  No.  3  (4.6  cm3). 
The  mean  signal  for  the  three  channels  was  calculated  and  is 
plotted  by  the  thicker  trace,  (b)  Correlation  between  Ap02  at 
individual  locations  in  the  tumors  in  response  to  carbogen  or  oxy¬ 
gen  for  the  five  tumors  in  Group  1  {R  >  0.8). 


4.  Discussion 

In  this  study,  we  have  simultaneously  measured  the 
arterial  Sa02,  the  global  changes  in  the  A[Hb02]  of 
tumor  vasculature,  and  the  regional  changes  in  the 
Ap02  of  tumor  tissue,  in  response  to  hyperoxic  (i.e., 
carbogen  and  oxygen)  gas  interventions  with  a  pulse 
oximeter,  an  NIRS  system  and  a  multichannel,  fiber¬ 
optic,  oxygen-sensing  system,  respectively.  All 
three  oxygen-sensitive  indicators  displayed  similar 


Fig.  7.  Correlation  between  mean  Ap02  and  A[Hb02]  for  the  five 
breast  tumors  ( R  >  0.86):  ♦,  transition  from  air  to  carbogen;  A, 
transition  from  air  to  oxygen. 


dynamic  tendency  in  response  to  carbogen  interven¬ 
tion  (Fig.  2). 

The  simultaneous  measurements  demonstrate  the 
compatibility  of  the  NIRS  system  with  the  FOXY 
fiber-optic  oxygen-sensing  system,  without  interfer¬ 
ence.  Both  systems  are  relatively  inexpensive  and 
provide  real-time  measurements,  but  the  multichan¬ 
nel  FOXY  fiber-optic  system  monitors  Ap02  in  spe¬ 
cific  locations,  whereas  the  NIRS  system  provides 
global  measurements.  Whether  AHb02  determined 
with  this  methodology  will  be  a  clinically  useful  pre¬ 
dictor  for  tumor  response  to  oxygen-dependent  inter¬ 
ventions  and  therapies  remains  to  be  determined. 
However,  it  is  established  that  measurements  of  p02 
have  prognostic  value  in  the  clinic18  20  so  that  corre¬ 
lations  between  p02  and  NIR  measurements  would 
be  very  important. 

We  have  previously  applied  a  polarographic  oxygen 
electrode  simultaneously  with  NIR.25  However, 
that  study  provided  only  a  single  local  p02  value,  and 
in  some  cases  correlations  with  global  NIR  measure¬ 
ments  were  very  poor.  The  optical  fiber  system  used 
here  allows  multiple  locations  to  be  interrogated  si¬ 
multaneously.  The  device  can  be  expanded  to  many 
channels,  but  our  system  uses  four  channels.  Un¬ 
fortunately,  probes  are  fragile,  and  the  oxygen- 
sensitive  coating  on  the  tips  is  readily  damaged. 
Thus,  we  only  had  three  probes  available  for  this 
study.  Indeed,  fiber-optic  probe  fragility  is  a  well- 
recognized  problem,  and  our  previous  experience 
with  the  more  expensive  OxyLite  system  was  also 
restricted  to  three  channels  owing  to  probe  damage.26 

The  FOXY  system  (~$13k)  is  much  less  expensive 
than  the  OxyLite  (~$48k),  and  its  mode  of  action  is 
also  simpler,  detecting  fluorescent  signal  intensity 
rather  than  fluorescence  lifetime.  It  seems  capable 
of  measuring  p02  across  the  whole  range  of  atmo¬ 
spheric  oxygen  tensions  (0-760  Torr),  whereas  the 
OxyLite  is  restricted  and  becomes  very  insensitive 
above  approximately  100  Torr.  However,  our  expe¬ 
rience  shows  that  although  the  FOXY  system  pro¬ 
vides  precise  measurements  of  Ap02,  absolute  values 
of  p02  may  not  be  reliable.  We  continue  to  perform 
validation  experiments.  By  contrast,  the  OxyLite 
system  seems  to  give  very  accurate  p02  values. 

Our  experience  shows  that  the  FOXY  probes  are 
much  easier  to  use  than  electrodes,  particularly,  in 
terms  of  calibration  and  stability.  Since  the  probes 
are  fragile,  we  insert  them  into  tumors  through  a  fine 
needle  (25  gauge),  which  readily  punctures  the  sur¬ 
rounding  skin  and  penetrates  tough  fibrous  tissues. 
The  needle  is  then  backed  up  from  the  tip  to  facilitate 
measurements.  The  probes  often  require  a  few  min¬ 
utes  to  settle  at  a  stable  baseline  value,  but  then  show 
good  baseline  stability  and  rapid  response  to  inter¬ 
ventions  [Figs.  2  and  6(a)].  They  are  easily  moved 
within  the  tumor  to  locate  regions,  presenting  a  par¬ 
ticular  p02  of  interest,  e.g.  hypoxic  or  well  oxygen¬ 
ated.  In  the  search  for  appropriate  locations,  probes 
are  moved  forward  to  interrogate  fresh  tissues  rather 
than  in  reverse,  since  blood  may  pool  in  the  tracks 
owing  to  vascular  damage.  However,  we  observe 
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minimal  bleeding  on  removal  of  the  probes  from  the 
tumors. 

We  have  found  no  interference  between  the  NIR 
and  FOXY  instruments,  although  any  tumor  motion 
associated  with  moving  the  fiber  probes  can  alter  the 
optical  contact  of  the  NIR  optrodes,  and  thus,  alter 
apparent  AHb02.  Thus,  baseline  AHb02  is  deter¬ 
mined  once  the  fiber  probes  are  situated.  New  fiber¬ 
optic  probes  of  the  FOXY  system  have  a  thick  coating 
of  fluorescent  gel  and  a  black  covering,  but  this  wears 
with  use  and  gradually  allows  reception  of  the  NIR 
light.  Since  the  LEDs  of  the  two  systems  operate  at 
very  different  wavelengths,  viz.  475  versus  760  nm, 
there  is  no  interference  for  detection.  The  detection 
of  local  NIR  light  by  the  FOXY  probe  opens  the  ex¬ 
citing  possibility  of  detecting  regional  hemoglobin  ox¬ 
ygen  saturation.  We  believe  the  FOXY  optical 
probes  could  be  moved  within  the  tumor  to  map  the 
distribution  and  path  of  the  transmitted  NIR  light, 
helping  to  explore  and  validate  the  optical  character¬ 
istics  of  the  tumor.  This  can  simultaneously  provide 
a  correlation  between  local  AHb02  and  Ap02. 

In  this  study,  we  have  examined  a  much  larger 
group  of  rats  than  previously.16*25  We  have  now 
shown  rigorously  that  the  two  hyperoxic  gases  induce 
similar  changes  in  vascular  oxygenation  (NIR)  and 
regional  tissue  p02  (FOXY)  in  this  type  of  rat  breast 
tumor.  These  data  are  consistent  with  our  previous 
observations  using  19F  NMR  imaging  (FREDOM)33  in 
this  tumor  type  and  also  in  rat  prostate  tumors.34*35 
If  the  two  gases  are  indeed  equivalent  in  terms  of 
manipulation  of  tumor  oxygenation,  it  could  have 
great  therapeutic  benefit  since  the  popular  carbogen, 
which  is  in  use  in  clinical  trials,36  can  cause  respira¬ 
tory  discomfort. 

The  current  data  show  that  AHb02  and  Ap02  are 
correlated  (Fig.  7),  and  thus,  such  noninvasive  obser¬ 
vations  could  have  value  in  the  clinic.  The  major 
deficiency  in  our  current  NIR  approach  is  lack  of 
spatial  discrimination,  and  thus  efforts  to  implement 
NIR  imaging  will  be  of  great  value.  It  will  also  be 
interesting  to  correlate  other  measurements,  such  as 
blood-oxygen-level-dependent  (BOLD)  proton  mag¬ 
netic  resonance  imaging,  which  provide  high  spatial 
resolution,  but  which  are  sensitive  to  vascular  flow 
and  volume  as  well  as  oxygenation.37 

The  biphasic  response  of  AHb02  to  carbogen  is  in¬ 
triguing,  and  we  believe  it  represents  the  distinct 
vascular  compartments  of  arterioles  (high  flow)  and 
capillaries.  However,  the  change  to  monophasic  be¬ 
havior,  when  carbogen  is  administered  second,  re¬ 
quires  further  exploration;  in  the  future,  we  propose 
to  test  various  concentrations  of  oxygen  and  carbon 
dioxide  and  air  to  separate  the  components  of  the 
response.  The  carbogen  dioxide  component  of  car¬ 
bogen  is  known  to  be  vasoactive;  however,  the  specific 
effects  may  depend  on  tumor  type,  site  of  growth,  and 
other  factors.9*38 

In  terms  of  vascular  oxygen  delivery,  the  data  in 
Table  1  reveal  the  progressive  movement  of  oxygen: 
td  (Sa02)  <  td  (A[Hb02])  <  td  (Ap02).  As  expected, 
switching  to  hyperoxic  gas  caused  the  systemic  arte¬ 


rial  Sa02  to  increase,  as  a  result  of  the  immediate 
combination  of  deoxyhemoglobin  with  oxygen.  The 
highly  oxygenated  blood  circulated  in  the  systemic 
vasculature  of  the  rats  (including  the  capillary  bed  of 
the  tumor  tissue),  resulting  in  a  delayed  increase  in 
[Hb02]  in  the  tumor  vasculature,  and  led  to  an  un¬ 
loading  of  oxygen  to  the  tumor  tissue.  For  the  biex¬ 
ponential  model  of  A[Hb02],  the  fast  component  has 
a  similar  time  constant  to  the  Sa02  measured  with 
the  pulse  oximeter  on  the  hind  leg,  strongly  suggest¬ 
ing  that  it  represents  arteriolar  oxygenation  in  the 
tumor.  In  this  study,  tumor  volumes  do  not  show 
any  direct  relation  with  time  constants  or  changes  of 
amplitude  in  response  to  hyperoxic  gas  interventions. 

It  is  increasingly  evident  that  oxygen  and  hypoxia 
play  important  roles  in  tumor  development  and  re¬ 
sponse  to  therapy.18  NIR  offers  an  attractive  non¬ 
invasive  means  of  investigating  tumor  oxygenation, 
particularly  in  terms  of  dynamic  response  to  inter¬ 
ventions,  but  we  had  previously  shown  a  potential 
mismatch  between  global  AHb02  and  local  Ap02.25 
The  data  presented  here  indicate  a  correlation  be¬ 
tween  the  global  NIR  measurements  and  mean  p02 
values  with  even  as  few  as  three  representative  loca¬ 
tions  per  tumor.  This  does  suggest  that  it  will  be 
important  to  develop  regional  NIR  measurements 
and  that  even  relatively  crude  mapping  could  reveal 
heterogeneity.  In  the  meantime,  we  believe  these 
studies  provide  further  evidence  for  the  value  of 
NIRS  to  explore  tumor  physiology. 

This  study  was  supported  in  part  by  the  Depart¬ 
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(VB).  We  are  grateful  to  Mengna  Xia  and  Dawen 
Zhao  for  their  assistance  with  data  processing.  We 
gratefully  acknowledge  Weina  Cui  for  helpful  discus¬ 
sions. 
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Appendix  5 


Investigation  of  breast  tumor  hemodynamics  using 
tumor  vascular  phantoms  and  FEM  simulations 
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Abstract:  We  have  developed  dynamic  tumor  vascular  phantoms  and  utilized  the  finite  element 
method  to  investigate  the  bi-phasic  behavior  of  oxyhemoglobin  concentration  increase  that  were 
observed  in  vivo  from  rat  breast  tumors  during  carbogen  inhalation. 

©2004  Optical  Society  of  America 

OCIS  codes:  (170.1470)  blood/tissue  constituent  monitoring;  (170.3660)  Light  propagation  in  tissues;  (170.4580)  Optical 
diagnostics  for  medicine;  (170.6510)  Spectroscopy,  tissue  diagnostics 

1.  Introduction 

Solid  tumors  develop  regions  of  hypoxia  during  their  growth  due  to  an  imbalance  between  the  rate  of  tumor  cell 
proliferation  and  branching  of  the  blood  vessels  [1,2].  Tumor  hypoxia  can  contribute  to  the  failure  of  radiotherapy 
[3],  some  forms  of  chemotherapy  [4]  and  photodynamic  therapy  [5].  In  addition,  a  number  of  clinical  trials  have 
found  that  patient  survival  measured  either  as  tumor  regression  or  as  local  control  depends  highly  on  tumor 
oxygenation  [6].  Therefore,  increasing  tumor  oxygenation  could  improve  cancer  therapy  efficiency.  As  one  means  to 
improve  tumor  oxygenation,  breathing  hyperoxic  gas  has  been  used  to  enhance  cancer  treatment  efficacy. 

Based  on  near  infrared  spectroscopy  (NIRS)  measurements,  our  previous  in  vivo  animal  study  has  clearly 
demonstrated  that  carbogen  (95%  C02  and  5%  02)  inhalation  can  improve  the  vascular  oxygen  level  of  breast 
tumors  [7].  The  observed  changes  in  tumor  vascular  oxygenation  (A[Hb02])  exhibited  a  rapid  increase,  followed  by 
a  gradual  but  significant  increase,  in  response  to  carbogen  intervention.  To  explain  these  biphasic  behaviors  of 
tumor  hemodynamics,  we  established  a  mathematical  model  based  on  Kety’s  approach  [8],  and  we  formed  a 
hypothesis  that  tumor  vasculature  is  comprised  with  a  well-perfused  and  poorly  perfused  region  to  explain  why  there 
are  two  different  time  constants  in  the  A[Hb02]  data.  The  mathematical  model  further  allows  us  to  associate  the 
signal  amplitudes  to  the  ratio  of  vascular  coefficients  and  the  ratio  of  the  perfusion  rates  in  the  two  different  regions. 

Even  though  we  associated  our  fitting  parameters  of  the  NIR  measurements  mathematically  with  physiological 
factors,  such  as  tumor  vascular  coefficients  and  perfusion  rates,  we  did  not  provide  solid  validation  for  the  model  at 
the  time.  In  this  study,  to  further  investigate  breast  tumor  vasculature  and  to  interpret  the  mathematical  model,  we 
designed  and  developed  a  tumor  hemodynamic  phantom  for  multi-channel  NIRS  experiments.  Furthermore,  we 
used  the  FEMLAB  software  to  simulate  the  dynamic  processes  of  light  distribution  in  tumor  hemodynamic 
phantoms,  giving  rise  to  better  understanding  of  the  association  between  the  NIR  readings  and  tumor  physiological 
parameters. 

2.  Materials  and  methods 

2.  /  Mathematical  Model  of  Tumor  Vascular  Oxygenation 

The  mathematical  model  for  tumor  vascular  oxygenation  during  carbogen  inhalation  has  been  described  in  details  in 
our  previous  report  [7].  Briefly,  we  followed  an  approach  used  to  measure  regional  cerebral  blood  flow  (rCBF)  with 
diffusible  radiotracers,  as  originally  developed  by  Kety  [8]  in  the  1950’s.  By  applying  Fick’s  principle  and  assuming 
that  tumor  has  two  distinct  perfusion  regions,  the  measured  NIRS  signal,  which  has  a  double-exponential  feature 
observed  in  the  rat  tumor  experiments,  can  be  expressed  as 

=  /,Ho[  1  -exp(-/i//  /,)]  +  KHo[l-exp(-/2//tt)] 

=  Aifl-expC-z/x,)]  +  A2[l-exp(-//t2)],  (1) 

where  Hc  is  the  arterial  oxygenation  input,/!  and  yx  are  the  blood  perfusion  rate  and  the  vasculature  coefficient  in 
region  1  for  the  well  perfused  region,  respectively;/  and  yi  have  the  same  respective  meanings  in  region  2  for  the 
poorly  perfused  region,  and  A{  =  ^Ho,  A2  =  #He,  Ti  =  yjf\,  t2  =  #//  Basically,  Au  A2,  Tj,  and  x2  can  be  determined 
from  the  NIRS  measurements  by  fitting  equation  (1)  with  the  experimental  data.  It  follows  that  the  ratios  of  two 
vasculature  coefficients,  yi/y2,  and  the  two  blood  perfusion  rates,/;//,  can  be  obtained  as 
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Zl=A  (2) 

Yl  a2  '  fl  r,/r2 

2.2  Tumor  Hemodynamic  Phantom  and  Multi-Channel  NJR  Spectroscopy 

To  represent  vascular  blood  vessels  in  tumor,  we  designed  a  vascular  modeling  device  (VMD)  by  winding  a 
small  diameter  tube  (I.D.=0.51  mm)  around  a  big  diameter  core  tube  (O.D.=14.4  mm).  Tumor  hemodynamic 
phantom  was  fabricated  by  embedding  two  VMDs  into  a  cylindrical  shape  of  gel,  which  is  made  of  Intralipid  and 
gelatin  powder  to  represent  tissues.  The  optical  properties  of  tumor  gel  phantom  were  measured  using  near  infrared 
tissue  spectrometer  (model:  96208,  ISS  Inc.,  Champaign,  IL),  and  those  values  were  close  to  tissue  optical  properties 
(pa=  0.032  cm'1,  ps’=9.2  cm'1  at  750  nm). 

Schematic  experimental  setup  for  tumor  hemodynamic  phantom  study  is  shown  in  Figure  la.  As  shown  here, 
NIR  light  at  730  nm  from  a  multi-channel  NIRS  system  was  delivered  to  the  tumor  hemodynamic  phantom,  and 
three  optical  detectors  were  placed  on  the  side  of  cylindrical  phantom  (Dl,  D2,  and  D3  in  Fig.  lb)  to  collect  signals, 
which  were  then  processed  in  the  computer.  A  certain  amount  of  ink  was  injected  into  the  two  VMDs  by  two  syringe 
pumps  (model  KDS200,  KdScientific  Inc.,  New  Hope,  PA),  and  the  ink  wastes  were  collected  in  the  waste  beaker. 
We  used  two  syringe  pumps  so  that  we  could  control  different  ink  flow  rates  for  different  VMDs. 


Fig.  1(a)  Schematic  diagram  of  tumor  vascular  phantom  experimental  setup,  (b)  The  geometry  of  FHM  simulation  for 
tumor  hemodynamic  phantom. 


2.3  FEM  Simulation 


To  further  investigate  bi-phasic  features  of  tumor  hemodynamics,  we  used  the  FEMLAB  software  (COMSOL 
Inc.  Burlington,  MA),  which  is  based  on  the  finite  element  method,  to  simulate  our  dynamic  tumor  phantom 
experiments.  The  geometry  of  FEM  simulations  for  the  tumor  hemodynamic  phantom  is  shown  also  in  Figure  lb, 
and  it  represents  a  2-dimensional  horizontal  cross  section  of  our  tumor  vascular  phantom.  The  large  circle 
represents  the  body  of  tumor  phantom  (diameter  =  4  cm),  and  each  VMD  used  in  the  dynamic  phantom  is 
represented  by  two  thin  rectangles  (0.1  x  2.4  cm  rectangle)  within  the  circle  as  a  blood  vessel  network:  left  two 
rectangles  are  for  the  vessels  with  a  fast  flow  rate,  and  right  two  rectangles  denote  the  vessels  with  a  slow  flow  rate. 

In  the  simulation,  we  used  a  uniform  reduced  scattering  coefficient  of  10  cm'1  for  both  the  background  and 
vasculature.  Two  absorption  coefficients  of  0.03  and  1.5  cm'1  were  used  to  simulate  the  phantom  and  ink  absorption 
in  VMDs,  respectively.  The  FEM  model  was  performed  with  1 147  elements  and  609  nodes  and  solved  for  the  light 
distribution  within  the  vascular  phantom  dynamically,  using  the  stationary  nonlinear  solver  type. 

3.  Results 


Figure  2a  shows  the  changes  of  absorption  from  the  three  detectors  on  the  tumor  hemodynamic  phantom.  Since  two 
VMDs  are  having  same  inner  diameter  (0.51  mm),  the  velocities  of  ink  in  VMDs  will  be  totally  depending  on  the 
ink  flow  rates  controlled  by  the  two  syringe  pumps.  The  first  step  for  this  experiment  was  injecting  ink  only  into  left 
VMD  with  a  flow  rate  of  20  ml/hr.  As  a  result,  D2  showed  largest  increase  of  Apa  while  Dl  showed  smallest 
increase  of  Apa  because  Dl  was  located  quite  away  from  the  left  VMD.  For  the  second  step,  right  VMD  was 
injected  with  ink  at  the  same  flow  rate  of  20  ml/hr.  Here,  Dl  showed  largest  increase  of  A(ia  and  D2  showed 
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smallest  increase  of  A|ja.  Then,  as  the  third  step,  the  ink  solution  was  injected  into  both  left  and  right  VMDs  with  the 
same  flow  rate  (20  ml/hr).  In  all  of  the  above  cases,  the  signals  from  D3  showed  similar  profiles  as  the  others, 
without  clear  bi-exponential  patterns.  For  the  fourth  step,  we  injected  ink  into  both  VMDs,  but  with  two  different 
flow  rates:  5  ml/hr  for  right  VMD  and  20  ml/hr  for  the  left  VMD.  Now  we  can  see  that  bi-exponential  behavior 
appears  on  D3  since  D3  is  detecting  signals  from  both  left  and  right  VMDs,  which  have  two  different  ink  flow  rates. 

For  the  further  analysis,  Apa  increases  detected  from  the  three  detectors  at  step  4,  as  shown  in  Figure  2a,  were 
fitted  by  using  our  mathematical  model  to  obtain  amplitudes  and  time  constants  values  (Figure  2b).  The  time 
constant  analysis  shows  that  the  Apa  increases  from  D1  and  D2  are  well  fitted  by  the  mono-exponential  model, 
while  Apa  increases  from  D3  is  fitted  well  with  the  double  exponential  model,  i.e.,  eq.  (1).  There  exist  a  fast  and 
slow  time  constants  with  Tj  =  18.1  ±  0.9  sec  and  x2  =  134.  ±  1 1  sec  from  D3,  and  they  are  well  matched  with  ij  in  D2 
(=  20.9  ±  0.5  sec)  and  n  in  D1  (=  131.  ±  3  sec).  All  of  these  suggest  that  the  bi-phasic  Apa  signals  obtained  from  D3 
with  a  fast  and  slow  component  indeed  result  from  the  dynamic  changes  in  both  left  and  right  VMDs.  We  have 
further  confirmed  the  same  results  from  the  FEM  simulations  on  tumor  hemodynamic  phantom  experiments. 


Fig.  2  (a)  Effects  of  overdose  KCI  injection  on  tumor  vascular  A[Hb02]  and  A[Hb]toui  for  a  breast  tumor  (12.7  cm'),  (b) 
relationship  between  the  tumor  oxygen  consumption  rate  and  tumor  volume. 

4.  Discussion  and  Conclusions 


We  have  developed  dynamic  tumor  phantom  models  to  demonstrate  our  hypothesis  on  tumor  hemodynamics 
during  carbogen  inhalation.  We  believe  that  the  bi-phasic  feature  of  tumor  blood  oxygenation  increases  during 
carbogen  inhalation  is  highly  associated  with  distinct  vascular  structures  of  the  tumors,  which  are  composed  of  well- 
perfused  and  poorly  perfused  region.  We  found  that  the  two  time  constants  in  tumor  hemodynamic  models  can 
result  from  different  blood  flow  velocities  or  anything  that  can  cause  changes  in  blood  flow  velocities,  such  as  blood 
vessel  diameter  and  length.  Also,  the  FEM  simulation  based  on  light  diffusion  equation  was  performed  to  simulate 
the  tumor  phantom  experiments,  and  the  results  were  well  matched  with  the  experiments,  further  supporting  our 
previous  hypothesis  and  the  mathematical  model. 
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